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Preface 



Near-Space Remote Sensing: Potential and Challenges deals with the role of near- 
space vehicles in supplying a gap between satellites and airplanes for microwave 
remote sensing applications from a top-level system description, with an aim for 
further research. 

Near-space is defined as the atmospheric region from about 20 to 100 km 
altitude above the Earth's surface. Very few sensors are currently operating in 
near-space, because the atmosphere is too thin to support flying for most aircrafts 
and yet too thick to sustain orbit for satellites. Nevertheless, potential benefits for 
vehicles operating in near-space may include possible persistent monitoring and 
high revisiting frequency (revisit the same site in a short time interval) that are 
critical to some specific radar and navigation applications, but not accessible for 
current satellites and airplanes. 

There is a region in near-space where the average wind is <10 m/s; hence, 
persistent coverage and high flying speed can be obtained for the vehicles oper- 
ating in this region. Moreover, near-space vehicles are relatively low cost when 
compared to satellites and airplanes. Additionally, as near-space is below iono- 
sphere, therefore, there are no ionospheric scintillations that will significantly 
degrade microwave communication and navigation performance which explains 
why near-space has received much attention in recent years and why several types 
of near-space vehicles are being studied, developed, or employed. 

Near-space vehicles offer the long-term persistence traditionally provided by 
satellites while providing the fast responsiveness of airplanes. Near-space thus 
offers an opportunity for developing new synthetic aperture radar (SAR) remote 
sensing techniques. Several potential applications such as passive surveillance, 
reconnaissance, high-resolution wide-swath remote sensing, and ground moving 
targets indication (GMTI) are described in this book. 

Near-space vehicle-borne SARs cannot replace spaceborne and airborne radars, 
but they could provide more efficiently remote sensing functionality than space- 
borne SARs. Although near-space vehicles have much smaller coverage area than 
satellites due to their lower altitude, they can still offer a regional coverage of 
hundreds of kilometers and provide cost-effective remote sensing services. 
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Near-space vehicle-borne SARs could also extend remote sensing services into 
areas with limited or no access to spaceborne and airborne SARs. Therefore, given 
their operational flexibility, near-space vehicle-borne SARs may supply a gap 
between spaceborne and airborne SARs which is the reason why we appeal to the 
system engineering community for more publication and support on the research 
and development of near-space vehicle-borne SARs. 

With great pleasure I acknowledge the many people who have influenced my 
thinking and contributed to my knowledge. I express my deepest gratitude to Profs. 
Qicong Peng and Jingye Cai at the University of Electronic Science and Tech- 
nology of China. They provided me with unprecedented freedom to spend my time 
on almost any topic that stimulated my curiosity. I also thank the support and 
encouragement of Prof. Xiaowen Li at the Institute of Remote Sensing Applica- 
tions, Chinese Academy of Sciences. 

I also wish to acknowledge the support provided by the Doctoral Program of 
Higher Education for New Teachers under contract 200806141101, Fundamental 
Research Funds for the Central Universities under contract ZYGX2010J001, and 
the Open Funds of the State Laboratory of Remote Sensing Science under contract 
OFSLRSS201011, the National Key Laboratory of Millimeterwave Technology 
under contract K200914 and the Key Laboratory of Ocean Circulation and Waves, 
Chinese Academy of Sciences under contract KLOCAW1004. 

As always the support of my wife Ke Yang is gratefully acknowledged for her 
constant and gentle encouragement. 

Finally, I thank Ms. Becky Zhao and Na Xu from Springer for their wonderful 
help in the preparation and publication of this manuscript. 
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Chapter 1 
Introduction 



Abstract Near-space, the region between controlled commercial airspace and low- 
earth orbit (LEO), offers several capabilities that are critical to emerging remote 
sensing applications, but not accessible to current satellites and airplanes. In this 
chapter, we explain what near-space is and why it should be exploited for remote 
sensing. 

Keywords Near-space • Space definition • Near-space environment • Satellite • 
Airplane • Remote sensing 



1.1 Background 

Near-space was an area out of reach of airborne aircraft, yet below the effective 
satellite orbit. However, in recent years new technology and interest have brought 
about a reexamining of this place between air and space [1]. 

1.1.1 What is Near-Space 

Near-space is defined as the atmospheric region from about 20 km altitude to 100 km 
altitude above the Earth's surface [2] , as shown in Fig. 1.1. Note that the lower limit is 
not determined from operational considerations, but from the international controlled 
airspace altitude [3]. 

Traditionally very few sensors are currently operating in near-space because the 
atmosphere is too thin to support flying for airplanes and yet too thick to sustain 
orbit for satellites. However, evolutionary advances in several technologies have a 
revolutionary advance in capability. One technology is the power supplies includ- 
ing thin, light-weight solar cells, small, efficient fuel cells, and high-energy-density 
batteries; the miniaturization of electronics and exponential increase in computing 
power, enabling extremely capable sensors in very small and light-weight packages; 
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Fig. 1.1 Near-space 
definition and its advantages 
while compared to space 
including geosynchronous 
orbit (GEO), middle earth 
orbit (MEO) and LEO, and 
airplane 
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and very light-weight, strong, and flexible materials that can resist degradation under 
strong ultraviolet illumination and relatively impermeable to helium or hydrogen [4] . 

Another two emerging technologies are high-altitude buoyant lift systems and 
plasma thrust technology [2]. As altitude increases, electrodynamic power transfers 
into the air which can be used for propulsion, cooling, and control. This allows us 
to call upon many creative electromagnetic circuits for our designs. Plasma technol- 
ogy supports the high-altitude and prompt global strike programs pursued by many 
institutes or laboratories. 

Near-space thus has become an area of exceptional interest in the last several years. 
Possible real-world uses in communication, radar, and navigation applications have 
come to light. 



1.1.2 Near-Space Environment 



There are no clouds, thunderstorms, or precipitation in near-space. Although the air 
density in near-space is very low, wind is still an important environment factor. Wind 
in near-space varies with altitude, time of year, and latitude, generally increasing with 
both latitude and altitude. In higher near-space there is no appreciable wind. Note that 
sustained winds will have an effect, but guests are more important to be considered 
when designing the near-space vehicles that have requirements to maneuver. 

Another important advantage is that near-space vehicles operate below ionosphere. 
Ionosphere is a region containing several distinct layers of charged particles sur- 
rounding the earth. Its main layers occur at about 70-90, 95-140, and 200-400 km 
above the earth. Ionosphere affects the electromagnetic signals that pass through it. 
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Depending on frequency and direction of propagation, some signals are slowed, atten- 
uated, or bent slightly while others can be completely absorbed or bent. Spacecraft 
thus needs to be designed to mitigate operations in this highly-charged environment 
while near-space vehicles can avoid this cost. Therefore, in near-space there are no 
ionospheric scintillations that will significantly degrade microwave communication 
and navigation performance which explains why near-space has received much atten- 
tion in recent years and why several types of near-space vehicles are being studied, 
developed, or employed [5, 6]. 



1.1.3 Why Near-Space Remote Sensing 

Current spaceborne and airborne synthetic aperture radars (SARs) are placing an 
important role in remote sensing applications; however, even as good as they envi- 
sioned or employed, it is impossible for our limited non-geosynchronous earth orbit 
satellites to provide a staring presence on a timescale of days, weeks, or months over 
a selected target or area of interest. Even if we can launch a satellite for particular 
application, it would only be in view for a short time. As an example, most LEO 
satellites have a specific target in view for less than lOmin at a time. In contrast, 
conventional airplanes cannot fly too high because there is no sufficient oxygen to 
allow the engines to operate. Generally, satellites operate in the orbits above 200 km 
and air-breathing airplanes cannot routinely operate much higher than 1 8 km. 

Thus we have two gaps. The first is a gap in capability of persistent remote sensing 
observation. The second is a gap where there is little sensors in the altitude between 
air and space. To overcome these disadvantages, some alternative platforms should 
be employed [7]. Moreover, the requirements of stealth and robust survivability 
in military applications also call for new radar platforms other than satellites and 
airplanes. Fortunately, these aims can be simultaneously obtained by using near- 
space platforms, at a fractional cost of the traditional platforms [8]. 

In fact, by placing SAR transmitters and/or receivers inside near-space vehicles, 
many functionalities that are currently performed with satellites or airplanes could be 
performed more efficiently than conventional airborne and spaceborne SARs. Some 
works on the use of near-space vehicle-borne sensors for communication applications 
have been reported [9]. Some other works have been reported on near-space vehicles 
in radar and navigation applications [3, 10-13]. 



1.2 Outline of the Chapter 

In this book, we discussed the role of near-space vehicles in supplying a gap between 
satellites and airplanes for remote sensing applications from a top-level system 
description, with an aim for further research. Several potential applications and 
corresponding challenges are described. The remaining chapters are organized as 
follows. 
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Chapter 2 investigates the advantages of near-space vehicles in remote sensing, 
while compared to satellites and airplanes. Near-space vehicles can be classified into 
free-floater, steered free-floater, and maneuvering vehicle. Each category is further 
explained in great detail. 

Chapter 3 is devoted to near-space vehicles in passive microwave remote sensing. 
Applications of this technique to regional passive remote sensing and challenges of 
synchronization and motion compensations are given as examples. 

Chapter 4 deals with near-space vehicles in wide-swath remote sensing. The con- 
straints on high-resolution and wide-swath in SAR remote sensing are derived. Sev- 
eral multi-antenna or multi-aperture techniques, such as multi-channel in azimuth 
and multi-channel in elevation, are proposed. 

Chapter 5 proposes a scheme of near-space vehicle-borne multiple-input and 
multiple-output (MIMO) SAR for ground moving targets indication (GMTI) appli- 
cations. This approach is more effective and robust than the conventional displaced 
phase center antenna (DPCA) SAR-based GMTI solutions. 

Chapter 6 provides a concise summary of the near-space remote sensing and 
discusses the realistic issues and future work. 
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Chapter 2 

Near-Space Vehicles: Remote Sensing 

Advantages 



Abstract Near-space provides a promise for future remote sensing applications. 
Instead of concentrating pay loads, in this chapter we give a brief overview of the basic 
types of near-space vehicles currently in use, in active development, or envisioned. 
Their advantages, limitations, and vulnerabilities for microwave remote sensing are 
investigated. 

Keywords Near-space • Remote sensing • Near-space vehicle • Free-floater • 
Maneuvering vehicle • Persistent coverage 



2.1 Near-Space Vehicles 

Some of the near-space vehicles already exist and a great deal of new near-space vehi- 
cles are currently in prototype [1,2]. Near-space vehicles can be classified into three 
major categories [3]: (1) free-floaters, (2) steered free-floaters, and (3) maneuvering 
vehicles. Free-floaters are like rudimentary rafts where the speed and direction of 
travel are completely determined by the direction of the current. Steered free-floaters 
are like sailboats where the current has a large effect on their motion but they can 
steer within that current. Maneuvering vehicles are akin to steamships: they can go 
where they want and stay there for as long as they like. 



2.1.1 Free-Floaters 

Free-floaters are essentially large balloons that float with the wind. They are normally 
manufactured in two types: zero-pressure and super-pressure. Zero-pressure free- 
floaters are similar to weather and recreational balloons, which have a venting system 
which ensures that the pressure inside the balloon is the same as the surrounding 
atmosphere. They are less vulnerable to puncture, since significant amount of the 
lifting gas must diffuse out before the list is lost. Imaging an inflated, light-weight 
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Fig. 2.1 Example 
free-floater. Image courtesy 
USA Air Force Space Battle 
Lab 




plastic garment bag floating on the wind; even if there are many holes in such a bag, 
it still can float in the air for a long time. In contrast, super- pressure free-floaters 
have a higher pressure inside than outside. They are constructed by strong materials, 
making them relatively insensitive to puncture damage. 

Once launched, free-floaters are at the mercy of the existing winds. Limited steer- 
ing is possible by variable ballasting, causing them to float at different altitudes 
to take advantage of different wind directions and speeds. However, they have no 
station-keep ability because no active steering or propulsion techniques are employed 
in these platforms. Free-floaters have already demonstrated commercial viability as 
communication platforms, e.g. Fig. 2.1. They can lift payloads of tens to thousands 
of pounds to over 30 km, depending on their volume [4]. 

The biggest drawback of most free-floaters is their payloads generally cannot 
be recovered. In recent years, by encasing the payloads in a high-performance 
autonomous glider, expensive or sensitive payloads can be recovered safely and 
reused. When the floater approaches the maximum range of the glider, the glider is 
cut loose from the floater. The payload then autonomously glides back from hundreds 
of kilometers away, staying aloft for several hours before landing safely on ground. 
A variety of such gliders are available today, ranging from extremely inexpensive 
plastic gliders with limited payload capability to much more complex and capable 
composite gliders. 



2.1.2 Steered Free-Floaters 

Steered free-floaters also drift on the wind, but they are able to exploit the wind much 
like sailing ships to maneuvering almost at will. Sailing requires the vehicle to be 
immersed in two media moving at different speeds. A large balloon at high altitude 
moves at a different speed through the air than a wing suspended below the balloon at 
a different altitude. The air around the wing is moving at a different speed than the air 
pushing the balloon. The entire platform is then steered when the differential wind 
between the two parts of the platform enables the wing to become aerodynamically 
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Fig. 2.2 Example steered 
free-floater. Image courtesy 
USA Air Force Space Battle 
Lab 



Fig. 2.3 Example 
maneuvering vehicles. 
Image courtesy USA Naval 
Research Laboratory 




effective. Figure 2.2 gives an example of steered free-floater of the USA Air Force 
Space Battle Lab. Steered free-floater technology has been commercially mature and 
military deployment is imminent. 

Steered free-floaters can be navigated with fairly high degree of precision, gener- 
ally going with the flow of the prevailing latitudinal winds and being able to speed 
up, slow down, and move perpendicular to those winds. With the limited steering, 
these balloons can stay on station for a short time. Their payloads could be more 
complex than those flown on basic free-floaters. 



2.1.3 Maneuvering Vehicles 



Maneuvering vehicles has a means of propulsion and a means of control. The propul- 
sion can rely on fossil fuel, nuclear, or solar energy. The control can be attained 
through both aerodynamic and aerostatic means. They can maneuver and thus fly to 
and station-keep over a specific area of interest. Maneuvering vehicles are the func- 
tional cross between satellites and airplanes. They can provide large footprint and 
long-mission durations that are commonly associated with satellites and responsive- 
ness that is commonly associated with unmanned aerial vehicles. Figure 2.3 gives 
one conceptual version of the USA navy's high altitude airborne relay and router. 
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Fig. 2.4 Example tethered 
balloon-like aerostat 




Maneuvering vehicles can bring large payloads (up to 1 ,000 tons) that are enough 
to be useful for remote sensing applications and can be recovered for repair and reuse. 
Persistent coverage is also possible for maneuvering vehicles. One tethered balloon- 
borne radar system (e.g., Fig. 2.4) has been operationally used in an air defense and 
drug enforcement network operated by the USA Air Force. This network uses two 
sizes of aerostats which carry radars and operate at an altitude of 4.57 km feet. 

More importantly, not constrained by orbital mechanics like satellites and high- 
fuel consumption like airplanes, they can move at a speed as fast as 1,000-1,500 m/s. 
Thus, maneuvering vehicles are potentially the most useful type for the remote sens- 
ing applications that require fast revisiting frequency and high-resolution imaging 
performance. Maneuvering vehicles have also been viewed as a possible substitute 
for satellites supporting communications and other missions ([5-19]). 



2.2 State-of-the-Art 

Recently many projects focusing on near-space vehicles have been funded. In this 
section, the major near-space projects and players are briefly discussed. We focus 
here on their major achievements, including trials and demonstrations. 



2.2.1 In North America 

Sky Station was a North American commercial project consisting of a solar-powered 
aerostatic high-altitude platform system planned by Sky Station International [20]. 
This platform could be maintained geo-stationary at altitudes of approximately 2 1 km 
in the stratosphere, over major metropolitan areas. The average platform dimension 
could be 200 m long and 60 m in diameter at its widest point. Sky Station International 
planned to deploy at least 250 platforms, one about 21km above every major city 
in the world for wireless communication services. An overview of the Sky Station 
project can be found in Ref. [21]. 

Another North American commercial project is the Pathfinder, Pathfinder plus, 
HELIOS, and Sky Tower managed by USA AeroVironment Corporation. The first 
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Pathfinder was designed and fabricated in the early 1980s to support a classified pro- 
gram. It concluded that the required technology had not reached a level where ultra 
long-duration flight under solar power could be achieved. In 1998, Pathfinder was 
modified by NASA into the longer-winged Pathfinder Plus for civil telecommunica- 
tions. This aircraft achieved the record of a flight altitude of 24 km. The following 
HELIOS was capable of continuous flight for up to 6 months or more at 21km [22]. 
The Sky Tower aims to market high-altitude platforms to deliver a wide range of 
applications, such as fixed broadband communications, narrowband and broadcast 
communications. The last of the key technology work was completed in 2005, and 
practical flight and ground operations were demonstrated in 2006 [23]. 

The primary ongoing North American military effort to develop near-space vehi- 
cles is the advanced concept technology demonstration (ACTD), which was ini- 
tiated in 2003 to design, build, and test a high-altitude aerostat prototype that is 
able to operate unmanned, maintain a geostationary position over 2 1 km for up to 
6 months, generate its own power, and carry a variety of payloads. In 2003 Lockheed 
Martin Naval Electronics and Surveillance Systems was selected by the U.S. Missile 
Defense Agency to perform the ACTD Phase II contract award. The $40,000,000 
contract through June 2004 calls for the design of a solar-powered high-altitude air- 
ship. The airship is planned to have a mission life of one month, operating at 20 km, 
while providing 10 kW of power to a 4,000 pound pay load. It is intended that it will 
become a part of the Ballistic Missile Defense System Test bed. Some of the early 
activities demonstrated the capabilities of station-keeping and autonomous flight. It 
will be used for military and civilian activities including: (1) weather and environ- 
ment monitoring, (2) short- and long-range missile warning, (3) surveillance, and 
(4) target acquisition. The ACTD program documentation summarizes the effort as 
having "some technical risk" but "enormous potential benefits". 



2.2.2 In Europe 

In Europe, mainly two organizations have funded research activities throughout the 
continent, the European Space Agency (ESA) and the European Commission (EC). 
The ESA was one of the first research institutions to promote research on near-space 
vehicles in the continent. The EC has funded projects to develop and make demon- 
strators for near-space vehicles and to promote strategic studies related to the future 
of this technology. The representative research projects include HeliNet (network of 
stratospheric platforms for traffic monitoring, environmental surveillance, and broad- 
band services), CAPANINA (broadband communications technology), UAVNET 
(unmanned air vehicles network), CAPECON (civil UAV applications and economic 
affectivity of potential configuration solutions), and USICO (UAV safety issues for 
civil operations). 

In 1998 the ESA awarded to the Lindstrand Balloons Ltd a design contract 
for a geostationary stratospheric unmanned airship. This airship should fly in 
the stratosphere at 21km altitude and could remain stationary over its intended 
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position [24] . The airship was designed to carry a 600 kg pay load with a relay station, 
surveillance radar, and a weather radar or sensor package. It mainly devoted to the 
provisions of civil services in synergy or competition with both satellite and terres- 
trial systems. The STRATOS was another ESA-funded research project. The main 
object was to explore the grounds for the development and operation of a European 
near-space vehicle system based on the performance of a conceptual design for the 
best suited stratospheric platform concept answering the needs of future telecommu- 
nication markets [25]. 

The HeliNet was a project based upon high altitude very long endurance unmanned 
solar aerodynamic platforms, funded by the 5th Framework Programme of the 
European Union Commission. HeliNet was a global project carried out by a transna- 
tional and multi-sectorial partnership of research departments at universities and 
companies from Italy, Spain, UK, Slovenia, Hungary, and Switzerland. The proto- 
type was based on the design of an unmanned solar-powered aircraft, named Heliplat. 
Heliplat was tailored for long-endurance operations at an altitude of 17 km, support- 
ing a pay load of 1 00 kg and offering a power of 800 W. Apart from the Heliplat, three 
prototype applications were examined as well, i.e., broadband telecommunications, 
remote sensing, and navigation/location [26, 27]. After the successful completion 
of the HeliNet project in 2003, the EC started the CAPANINA research project, 
which is being partially funded by the 6th European Union's Framework initia- 
tive. Built on the HeliNet project, CAPANINA aims at the development of low-cost 
broadband technology from near-space vehicles to deliver cost-effective solutions to 
users in urban and remote rural areas, or to users traveling inside high-speed pub- 
lic transport vehicles at speeds up to 300km/h. One of the greatest achievements of 
CAPANINA was the fact that it proved credibility in the use of near-space vehicles to 
deliver broadband services. Three trials were completed successfully. The first trial 
was conducted in UK in 2004 using a very low-altitude tethered airship platform at 
300 m altitude. The second trial was conducted and completed in Sweden in 2005 
using a free-floater that could reach 25 km altitude. The last trial was conducted in 
USA in 2007 using a free-flight stratospheric balloon. 



2.2.3 In Asia-Pacific 

Asia-Pacific projects and activities on near-space vehicles were mainly undertaken by 
Japan, Korea, China, Australia, New Zealand, and Malaysia. Skynet was one Japanese 
project for the development of a balloon based on stratospheric platforms capable 
of operating at an altitude of 20 km and carrying on-board mission payloads for 
communications, broadcasting, and environmental observations. In the framework 
of this project, two prototype airships have been developed. The first airship, named 
Ground-to-stratosphere, has no propulsion system, and was successfully used to 
obtain thermal, buoyancy, and position control. It can ascend to the altitude of 15 km 
and descend to a planned area in the ocean. The second airship, named low-altitude- 
stationary, has a propulsion system. It will be used in order to obtain station-keeping 
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mechanisms. The propulsive propellers were mounted on both the stem and stern of 
the airship. A solar power subsystem of solar cells and regenerative fuel cells was 
provided to supply a day-night cycle of electricity for airship propulsion [28]. 

In Korea research activities on near-space vehicles were conducted mainly by 
ETRI (Electronics and Telecommunications Research Institute), whereas KARI 
(Korean Aerospace Research Institute) deals with airship research and development 
[29]. The main aim is to develop an unmanned stratospheric airship and ground sys- 
tems for basic operation and control of the airship. This research program consists 
of three phases. In the first phase one 50 m scale-size unmanned airship was built. 
The second phase aims to demonstrate the feasibility of developing a stratospheric 
platform that could be used at around 20 km altitude. The third phase aims to develop 
a full-scale 200 m airship that could carry telecommunications and remote sensing 
payloads weight up to 1,000 kg. 

In 2007 international agreements were achieved among the Malaysian govern- 
ment, QucomHaps Co. Ireland and the proprietor of Russian M-55GN stratospheric 
aircraft. The main aim was to provide cost-effective nation-wide wireless access 
to broadband connectivity using M-55GN stratospheric aircraft. The M-55GN is a 
piloted plane, which can fly in a circular corridor at an altitude of approximately 
21 km with a flight endurance of approximately 5 h. It is an all-weather single-seater 
stratospheric aircraft capable of operating both day and night, even in critical environ- 
mental conditions and strong cross-winds. Additionally, Chinese National Natural 
Science Foundation have also funded several research projects on near-space vehicles 
[30-33]. 



2.3 Comparative Advantages 

While compared to current satellite and airplane platforms, near-space vehicle plat- 
forms have many superiorities for microwave remote sensing applications [34]: 



2.3.1 Inherent Survivability 

Near-space vehicles, especially free-floaters, have an inherent survivability. Free- 
floaters have extremely small radar cross and thermal cross-sections, making them 
relatively invulnerable to most traditional tracking and positioning methods. Esti- 
mates of free-floater's radar cross-section (RCS) are on the order of hundredths of a 
square meter. In fact, at near-space altitude, free-floaters will be small optical targets, 
only showing up well when the background is much darker than them. Consequently, 
the acquisition and tracking problem will be very difficult, even without considering 
what sort of weapon could reach them. Surface-to-air missiles (SAMs) may be a 
threat, but they are most likely not designed to engage a non-maneuvering target at 
that altitude. Economics also discourage such an exchange because the trade between 
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Table 2.1 The maximum observation time (min:s) for selected LEO satellites at different orbital 
altitude (km) and different angle (degrees) above horizon 
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200 7:49 
300 9:35 
400 11:10 




5:37 
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2:24 1:27 
3:08 1:54 



an inexpensive near-space vehicle and a missile is cost-prohibitive. In modern elec- 
tronic battlefield, enemies may only want to launch their missiles for higher value 
and immediately destroyable targets. Free-floaters floating across their territory in 
numbers may not be worth of the effort and expense to shoot down. Moreover, even if 
the acquisition and tracking problems are overcome, near-space vehicles are difficult 
to destroy. 

On the other hand, there are defensive options available to help near-space vehi- 
cle's survivability. One such option could be deceptive "chaff" similar to that carried 
on modern fighter aircraft. This chaff could be dispensed in hopes of confusing radar 
guidance. Also, using vehicle decoys is an option. The vehicles themselves are rel- 
atively inexpensive and a simple fake payload could be attached. They might have 
imitated electronic and infrared signatures, further making it difficult to discriminate 
between the vehicles. 



2.3.2 Persistent Region Coverage or High-Revisiting Frequency 

The most useful and unique aspect of near-space vehicles is their ability to provide 
persistent region coverage or high-revisiting frequency. Space technologies have sig- 
nificantly revolutionized modern battlefield and remote sensing; however, persistent 
coverage, which is highly desired, is still unavailable through satellites or airplanes. 
Table 2.1 shows the observation time for selected LEO orbits. Air-breathing air- 
planes provide responsive persistence for the duration of their limited loiter times. 
The longest persistence that we can currently expect from an air-breathing asset is 
about a day or so. 

Fortunately, persistent coverage can be achieved by using near-space vehicles. 
The altitude of near-space is above troposphere and atmosphere region where most 
weather occurs. There is no cloud, thunderstorm, or precipitation. Moreover, propul- 
sion technique can be further applied to countermine possible mild winds in near- 
space. Near-space free-floaters can thus stay at a specific site for a long time, a 
persistence is thus possible. This advantage has a particular value for the applica- 
tions that require persistent monitoring. 

Near-space maneuvering vehicles can use a variety of schemes for propul- 
sion, including conventional propellers and unconventional buoyancy-modification 
schemes that allow the vehicles to propel themselves. Maneuvering vehicles are 
the functional cross between satellites and airplanes. They can fast-fly or station- 
keep over a specific position, to provide large footprint and long-mission duration. 
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Fig. 2.5 Ground coverage 
area as a function of 
looking-down angle for 
different flying altitude 
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The near-space vehicles currently receiving technology demonstration funding will 
be able to stay on station for 6 months, and planned follows are projected to stay aloft 
for years [4]. 

2.3.3 Relative High Sensitivity and Large Footprint 

Near-space vehicles are much closer to the targets than their orbital cousins. Distance 
is critical to receiving low-power signals. From the radar equation [35] we know that 
the received signal power attenuates at a square of the distance from the transmitter to 
the target, while that of an active antenna attenuates at a fourth power of the transmitter 
distance. Considering a point at nadir, near-space vehicles are 10-20 times closer to 
their targets than a typical 400 km LEO satellite. This distance differential implies 
that it could detect much weaker signals (10-13 dB weaker). 

On the other hand, near-space vehicles will have an impressive footprint from such 
high altitudes. Figure 2.5 shows the ground coverage area as a function of looking- 
down angle for different flying altitude. Although orbiting satellites traditionally have 
a much larger footprint than near-space vehicles, they do so at the cost of persistence 
and signal strength. Additionally, being lower than satellites also brings about another 
advantage to near-space vehicles: they fly below the ionosphere. Ionospheric scin- 
tillation is very difficult to predict, but it can disrupt satellite communication and 
navigation performance significantly. Fortunately, ionospheric scintillation has no 
impact on near-space vehicles. 

2.3.4 Low Cost 



The inherent simplicity, recoverability, and relative less requirements of complex 
infrastructure all contribute to this advantage. The cost of a free-floater and the 
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required helium to lift 100 pounds to an altitude of 20.42km is approximately $1,400. 
The cost of a larger free-floater and helium to lift 1 ,000 pounds up to the same altitude 
is approximately $10,000. The price of a high-altitude vehicle will be on the order 
of millions of dollars. The final goal of the Lockheed-Martin is to keep the cost per 
vehicle at $50-60 million. These costs can be compared to current UAV's such as the 
Predator ($4.5 million) and the Global Hawk ($48 million). Also, satellite expenses 
are traditionally enormous, with typical satellites costing roughly $60-300 million 
per copy. In addition to this cost is the expense to launch them into orbit, which 
adds another $10-40 million. Additionally, not being exposed to electronic radiation 
common to space environment, payloads flown in near-space vehicles require no 
costly space-hardening manufacture. 



2.4 Limitations and Vulnerabilities 

Near-space vehicles also have several limitations and vulnerabilities. The most 
prominent weaknesses are launch constraints, survivability constraints, and legal 
constraints. 



2.4.1 Launch Constraints 

Weather will be a risk factor that could be significant if vehicles are not furnished with 
reliable sensors for on-site meteorological data with which vehicle controllers can 
predict turbulence, icing, and violent gusts that jeopardize the craft. The experience 
with high-altitude tropospheric operation from around-the-world balloonist teams 
and weather teams must be collected and codified to aid computer predictions at 
higher altitudes. As a vehicle will be in the troposphere for over five hours while 
descending to its home-mooring base, the weather conditions will have to be within 
allowable parameters before a letdown can commence. This requirement could cause 
the vehicle to hold at near-space for up to two-five days before descending. Launch 
operations could cause similar delays. 

Note that satellites face similar launch constraints, but those constraints only have 
to be met during launch. UAVs and manned aircraft are also subject to similar launch 
and recovery constraints, although their limitations are less stringent than those for 
near-space vehicles. 



2.4.2 Survivability Constraints 

It is a hasty assumption to think that near-space vehicles will continue to remain out 
of harm's way in the future. As near-space technology develops, it should be assumed 
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that enemy technology will also develop, and air defense systems will improve. Near- 
space was not historically a viable region to operate and missiles were not developed 
to reach that high. However, a modern surface-to-air missile (e.g., SA-2) can reach 
as high as jet aircraft flew. The shoot-down of the U2 aircraft carrying Francis Gary 
Powers is an example showing that high-altitude flight is not impervious to surface 
weapons. 

An example of a fairly modern surface-to-air weapon is the the SA-10 Grumble. 
It is available for sale on the open market. Its newer versions boast a range of 56 
miles and an effective intercept ceiling of 26.82 km. This height extends well into 
that is recognized as near-space and exceeds the altitude of some of the vehicles 
forecast to operate there. Thus, the lower reaches of near-space are not beyond 
current conventional weapons. Near-space vehicles need to be flexible, replaceable, 
and cost-effective, otherwise their use may be very limited in military applications. 



2.4.3 Legal Constraints 

Freedom of overflight is another limitation. The legal status of the near-space regime 
is a gray area that is not directly addressed by treaty or policy. Near-space is not a 
new legal regime; the question is only whether it falls under air law, where nations 
claim sovereignty, or space law, where overflight rights exist. Due to lack of clear 
legal precedent governing the near-space regime, there is considerable disagreement 
over whether overflight rights exist [4]. 



2.5 Concluding Remarks 

Near-space has become an area of exceptional interest for space professionals in the 
last several years. Near-space vehicles can function as surrogate satellites but offer 
the advantage of shorter transmission distances for relaying ground-based commu- 
nication and ranges shorter than those of space-borne sensors for surveillance appli- 
cations. In military, the persistent surveillance from a fixed position by near-space 
vehicles, in contrast to periodic snapshots from the moving platforms that satellites 
or airplanes provide, allows continuous collection and comparison analysis over 
time of terrain covered by different sensors, such as infrared, electro-optical, and 
hyper-spectral imagery. In commercial, the comparative advantage of near-space 
vehicles to other means of broadband services has been assessed by many authors 
[1, 2, 7]. Commercial manufacturers are proposing near-space vehicles, including 
fixed-wing aircraft and high-altitude aircraft, to serve as surrogate satellites at a pre- 
sumably reduced cost. These efforts are driven by a desire to expand commercial 
high-bandwidth data services. The high cost of using satellites for that purpose has 
motivated the development of near-space vehicles for these commercial purposes. 
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Near-space vehicles can thus form an additional layer of persistence between 
satellites and airplanes. They improve upon the long-term persistence traditionally 
provided by satellites while providing the fast responsiveness of airplanes. Near- 
space vehicles thus offer an opportunity for developing new radar remote sensing 
techniques. First, they can support uniquely effective and economical operations. 
Second, they enable new remote sensing techniques. Third, they provide a crucial 
corridor for prompt regional strike. The details will be discussed in the following 
chapters. 
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Chapter 3 

Near- Space Vehicles in Passive 

Remote Sensing 



Abstract There is a region in near-space where the average wind is less than 10 m/s; 
hence, persistent coverage and high flying speed can be obtained for vehicles oper- 
ating in this region. For this reason, near-space has received much attention in the 
recent years. In this chapter, we consider mainly the role of near-space vehicles in 
passive remote sensing applications from a top-level system description. 

Keywords Near-space • Passive remote sensing • Surveillance and reconnaissance • 
Synthetic aperture radar (SAR) • Bistatic SAR • Persistently monitoring 



3.1 Near-Space Vehicles in Passive Surveillance 

Spaceborne and airborne radars have played an important role in surveillance and 
reconnaissance; however, regardless of how well they are envisioned and employed 
[1], it is impossible for limited spaceborne radars to provide a persistent coverage for 
an area of interest, because generally speaking, most LEO satellites have a specific 
target in view for less than ten minutes at a time and revisit the same site infrequently. 
Similarly, persistent coverage is also impossible for airborne radars. By contrast, per- 
sistent coverage is possible for near-space vehicles. Near-space vehicles, especially 
free-floaters, are inherently survivable [2], Moreover, several types of such vehicles 
are being studied, developed, or employed [3-6]. Their advantages provide a poten- 
tial to passive surveillance applications. However, little work on near-space vehicles 
in radar and navigation applications has been reported [7-10]. 



3.1.1 System Configuration 

The simplest and most natural application of near-space vehicles for surveillance 
is the passive radar system [11]. This system involves placing a passive receiver 
inside near-space vehicles and utilizing opportunistic illuminators such as global 
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Fig. 3.1 System 
configuration with 
two-channel receiver 
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positioning system (GPS) receivers and spaceborne or airborne imaging radars. The 
possibility of using GPS signals reflected off the sea surface and received by air- 
or spaceborne sensors have been proved by several authors [12-15]. This passive 
remote sensing can provide persistent monitoring without significantly impacting 
the normal day-to-day human activities [16]. 

As shown in Fig. 3.1, the passive receiver contains two channels: (1) one channel 
is used to receive the scattered signals with which detecting of targets is attempted; 
and (2) the other channel is employed for receiving the direct-path signals, which are 
used as the reference signal for subsequent matched filtering. This system offers two 
other advantages [17-19]: the first advantage is the potential of bistatic observation 
and the second is regional persistence. 

Another configuration using two or more passive receivers is also implementable. 
In this case, each receiver performs its own matched filtering. The total results can 
then be combined in some manner to provide a consistent detection or location. 
Additionally, a configuration using a single-channel receiver is also feasible. In this 
case, the received signals would contain the energy from both the direct-path channel 
and the scattered channel. Once they are separated, matched filtering can then be 
obtained successfully. 



3.1.2 Signal Models 



As an example, we consider the GPS LI signal [20] 



s L i(t) = <j2Pid(t)c(t) cos (2nf Ll t + <9 L1 ) + J2P^d(t)p(t) sm(2jr f Li t + B L[ ) 

(3.1) 
where P/ and Pq are the respective carrier power for the in-phase and quadrature- 
phase components, d(t) is the 50 bps data modulation, c{t) and p{t) are the respective 
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Fig. 3.2 Functional blocks 
of the correlation processor 
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C/A and P pseudorandom code waveforms, fn is the L\ carrier frequency in Hz 
per second, and On is a common phase shift in radians. The quadrature carrier Pq 
is approximately 3 dB less than Pj . Then, the received GPS scattered signal is 



X r (t) = a r SLl(f - X r ) 



(3.2) 



where a r and r r are the attenuation of the scattered signal and the corresponding 
delays in time, respectively. 

If matched filtering is applied with a reference function of the direct-path signal, 
the P codes will be masked by the C/A codes at the matched filtering output. This 
problem can be resolved by filtering out the C/A components in the direct-path chan- 
nel. The corresponding process is shown in Fig. 3.2. However, only a very limited 
resolution can be obtained from this method. Practically, GPS provides three types of 
measurements: pseudorange, carrier phase, and Doppler. Pseudorange measurement 
is based on the correlation of a satellite's transmitted code and the local receiver's 
reference code, which has not been corrected for errors in synchronization between 
the transmitter's clock and the receiver's clock; hence, it is a time error biased mea- 
surement. Carrier phase measurement is the difference between the phases of the 
receiver generated carrier signal and the carrier signal received from a satellite at the 
measurement instant. Carrier phase gives more precise measurements than pseudo- 
ranges, by estimating its instantaneous rate, or accumulated phase. 

Taking carrier phase measurement as an example, if the phase of the received 
carrier signal for the GPS satellite is denoted as (pg and the phase of the reference 
carrier signal generated by the receiver as <j>r , then there are 



<Ps(t s ) = 4>s(T s ) + </>s(Ar s ), <p R (t R ) = <Pr(T r ) + (p R (At R ) 



(3.3) 



The time is an epoch considered from an initial epoch. The corresponding carrier 
beat phase is 



4>%(f R ) 



4>s(t s ) ~ Mr) = 4>s(T s ) - <Pr(T r ) + f Ll At s - f Ll At R (3.4) 
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Fig. 3.3 Geometry model for target location 



The term fn/co with co the speed of light can be used to convert the geometric 
distance p into cycles. By accounting for all error sources affecting the carrier phase 
measurements, the observation in cycles can be represented by 



<Pm - P + Sp + c (At s - At R ) - d i( 



£0 



XN 



(3.5) 



where 8p is the orbital error, d{ 0n is the ionospheric advance, e^ is the receiver 
phase noise, X is the wavelength of the GPS carrier, and N is the number of integer 
cycles. The details can be found in many GPS textbooks. In this way, the transmitter- 
to-receiver distance and transmitter-target-receiver distance can be determined in a 
high-precision manner. 

3.1.3 Target Location 



To determine the instantaneous target position, we consider Fig. 3.3, which shows 
the passive bistatic radar geometry relations for target location. There are 



sin(ZTDA) 



R RC = L cos (90° -6 T ) 
R DF = R R cos (90° - 6 R ) 

h t -h r + R R sin(90° 
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where h, and h r are the transmitter altitude and the receiver altitude, respectively. 
The other variables are illustrated in Fig. 3.3. We can get 

r>2 j 2 p2 

K-r — U — An 

cos(g)= T * (3.9) 

ZLKr 

L 2 sm 2 (9 T ) + R 2 R sm 2 (9 R )-R 2 T cos 2 (ZTDA) 

cos(y) = . (3.10) 

2LR R sm(0 T ) &in(6 R ) 

Since there are L ^ Rr (target-to-receiver distance) andR j (transmitter-to-target 
distance) ^ Rr , we have 

cos (ZTDA)^^ (3.11) 

where the y and 6 are illustrated in Fig. 3.3. 
Suppose 

AR = R T + Rr- L (3.12) 

we can get 

_ { /AR(2L + AR) - 2R R (L + AR)\ 



cos 



cos 1 



2LR R ) 



2LR R sm(8 T )sm(0 R ) 



(, _ ( H-R R COS(8 R ) \ 2 \ 
y \ L+AR-R R ) J 



1/2 



(3.13) 

As the AR and L can be determined with GPS carrier phase measurement and 
correlation processor, the target-to-receiver distance Rr can be calculated from this 
equation. Next, the instantaneous target position (x s , y s , z s ) can be estimated by 

(x s , y s , z s ) = (x r - Rr cos (90° - Or) cos (e), 

y r - R R cos (90° - R ) sin(e), z r - Rr sin(90° - Or)) (3.14) 

where (x r , y r ,z r ) is the location of the receiver. 



3.1.4 Power Budget Analysis 

Using GPS satellite as an illuminator for near-space vehicle-borne passive bistatic 
radar presents a problem of signal detectability, since the reflected GPS signal is very 
weak. After coherent bistatic data processing, the final signal-to-noise ratio (SNR) 
can be represented by [21] 



24 



3 Near-Space Vehicles in Passive Remote Sensing 



Fig. 3.4 Example NESZ 
performance 
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(3.15) 



where Flo & 3 x 10 _14 Wt/m 2 , ero is the radar cross-section (RCS) parameter, K 
is the Boltzmann constant, 7b is the system noise temperature, F n is the noise figure, 
A r is the effective receiver antenna area, and r\ is the loss factor. 

A quantity directly related to radar detection performance is the noise equivalent 
sigma zero (NESZ), which is the mean RCS necessary to produce an S7V7? ra( i ar of 
unity. The NESZ can be interpreted as the smallest target RCS, which is detectable 
by the radar system against thermal noise. Setting SNR rd( i ai: — 1, Eq. 3.15 gives 



NESZ 



ATTRf.F„KT 
n A r T s r) 



(3.16) 



As an example, assuming a typical system with the following parameters: T s — 
1,000 s, 7b = 300, r\ = 0.5, and F„ = 2dB, the calculated NESZ is illustrated in 
Fig. 3.4. This result shows that an RCS requirement comparable to the current radar 
systems is possible for near-space vehicle-borne passive bistatic radar systems. 

3.2 Near-Space Vehicles in Passive BiSAR Imaging 



The use of the GPS illuminator has the advantages of entire planet coverage and 
simple transmitter-receiver synchronization [22], although high-resolution imaging 
cannot be obtained due to its limited signal bandwidth [23, 24] . In this section, we fur- 
ther discuss near-space vehicles in passive bistatic synthetic aperture radar (BiSAR) 
imaging when spaceborne or airborne radar is used as the opportunistic transmit- 
ter. Although the passive near-space vehicle-borne receiver may be stationary, an 
aperture synthesis can still be obtained by the transmitter's motion. 
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Fig. 3.5 General geometry 
of near-space vehicle-borne 
passive BiSAR 




BiSAR has received much recognition in the recent years [25]. But the published 
BiSAR systems are mainly azimuth-invariant configurations [26], for example the 
" Tandem" configuration where the transmitter and receiver move one after another 
with the same trajectory and the "translational invariant" configuration where the 
transmitter and receiver move along parallel trajectories with the same velocity. 
As near-space vehicle-borne receivers and spaceborne or airborne transmitters may 
have unequal velocities and/or unparallel trajectories, here, we consider mainly the 
azimuth-variant processing challenges and potential solutions. 



3.2.1 System Imaging Performance 

We consider one general near-space vehicle-borne passive BiSAR imaging geome- 
try, as shown in Fig. 3.5. The transmitter and passive receiver fly on a straight but 
nonparallel path with a constant but unequal velocity. The azimuth time is chosen to 
be zero at the composite beam center crossing time of the reference target. 



3.2.1.1 Imaging Time and Imaging Coverage 



For BiSAR imaging, both the transmit and receive antennas are steered to obtain 
an overlapping beam on the ground, and the exposure of a target is governed by 
this composite beam pattern. Thus, it is necessary to analyze the imaging time and 
imaging coverage. A space-time diagram-based approach was investigated in [27], 
but a rectangular approximation is employed. In fact, the whole imaging time of the 
near-space vehicle-borne BiSAR is determined by 



Dq,t + Da,r 

\v, -v r \ 



(3.17) 
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where v t and v r are the transmitter velocity and receiver velocity, D at and D ar 
are the illuminated ground coverage in azimuth direction for the transmitter and the 
receiver, respectively. They are determined, respectively, by 

h, /0 a ,\ h r (Bar\ 

D at ^2 — tan( — ), D a r « 2 — tan( — ) (3.18) 

sm(n) V 2 / sin(y,-) \ 2 / 

wherey f andy r are the incidence angle, and 9 a ,t and0 a r are the antenna beam width in 
azimuth direction for the transmitter and receiver, respectively. The imaging coverage 
in the azimuth direction can then be represented by 

Laz - D a.r ~ V r ■ T image (3.19) 

Similarly, the imaging coverage in the range direction is determined by 



<¥) 



L ra = D r ^2— %— t 
sin(y r ) 



where D rr and rj are the imaging ground coverage and antenna beamwidth in 
range direction for the receiver, respectively. 



3.2.1.2 Range Resolution 

The instantaneous range history of the transmitter and receiver to an arbitrary point 
target (x, y, 0) is 



R = J(x- x t ) 2 + (y- y t ) 2 + hj + J (x - x r ) 2 + (y - y r ) 2 + hj (3.21) 

where (x t , y t , h t ) and (x r , y r , h r ) are the coordinates of the transmitter and the 
receiver, respectively. We then have 

dR dR 

dx dy 

= [ sin(a r ) cos (£ f ) + sin(a r ) cos (f,-)]i x + [ sin(£,) + sin(^ r )]i y (3.22) 

where a t — oc t (x) and a r — a r (x) are the instantaneous looking-down angles, £ r = 
£t(x, y, yt) and C, r — C, r (x,y; y r ) (y t , y r is the instantaneous location in y-direction) 
are the instantaneous squint angles. There is 

\VR\ = V[sin(a,)cos(C,)+ sin(or r ) cos (£,.)] 2 + [sin(?,) + sin(^)] 2 (3.23) 

As the range resolution of a monostatic SAR is co/2B r with B r the transmitted 
signal bandwidth, the range resolution (in x direction) can then be derived as 

C °/ Br l n IAS 

Pr = ■ (3.24) 

[Vi?| frinfe,) 
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with 

sin(f f ) + sin(f r ) 



t; xy — arctan 



\ sin(i 



a, ) cos (£ f ) + sin(a,.) cos (£ r ) y 



Thus, we have 



[ sin(a,) cos (£,) + sin(a,-) cos (£>)] 
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(3.25) 



(3.26) 



We can conclude that the range resolution is determined by not only the transmitted 
signal bandwidth, but also the specific BiSAR configuration geometry. 



3.2.1.3 Azimuth Resolution 

At an azimuth time t, the range sum to an arbitrary reference point is 



R(r) = JR% + (v f r) 2 + J R* + (v r r) 



vo 



(3.27) 



where R,q and R,q are the closest ranges to a given point target when the transmitter 
and receiver move along their trajectories, respectively. The Doppler chirp rate is 
derived as 



kd(-c) = 



1 d 2 R(x) 
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Rto 



cos 



X 3 2 r 
The synthetic aperture time is determined by [28] 

XR,Q XR r Q 



( Vtt )+ V ^cos( V ^)] (3.28) 
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L t v t ' L,-v r 



(3.29) 



where L t and L r are the transmit antenna length and receive antenna length, 
respectively. 

— — ^ — — the corresponding Doppler bandwidth is 



If -j^ > -j 
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The azimuth resolution can then be expressed as 
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(3.31) 
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Table 3.1 Example BiSAR configuration parameters using near-space vehicle-borne receiver 
Parameters Transmitter A Receiver A Transmitter B Receiver B 



Carry frequency (GHz) 


5.33 


5.33 


1.25 


1.25 


Flying altitude (km) 


800 


20 


10 


20 


Flying velocity (m/s) 


7,450 


5 


100 


5 


Signal bandwidth (MHz) 


16 


16 


500 


500 


Beam incidence angle (°) 


30 


60 


60 


60 


Beamwidth in range (°) 


2.3 


15 


15 


15 


Beamwidth in azimuth (°) 


0.28 


15 


13.75 


13.75 


Imaging time (s) 




1.8670 




85.3248 


Imaging coverage (x,y)(km) 




(6.0808,6.0714) 




(2.3448, 3.0404) 



where 



' eq 



VJ + V? - 2v t V r COS (jT - £') (3.32) 



is the equivalent velocity, /J' is defined as the angle between the transmitter velocity 
vector and the receiver velocity vector. Similarly, if j^y- < ^— ^ , there is 



veq 
Pa = 



^sin 



(fe) + ^ in (^) 



(3.33) 



From Eqs. 3.31 and 3.33 we can see that, unlike the azimuth resolution of general 
monostatic SAR determined only by the azimuth antenna length, here the azimuth 
resolution will be impacted by its BiSAR configuration geometry and velocity. When 
R t0 — R r0 , v t — v r and L t — L r , there is p a — L t /2 — L r /2, In this case, it is just 
a monostatic SAR. If v t — or v r — and L t — L r , the azimuth resolution is found 
to be p a — L t — L r . This case is just a fixed-transmitter or fixed-receiver BiSAR 
[29]. Note that there is a local minimum for the azimuth resolution, which is L r /2 
or L,/2 (depending on whether ^ > ^ or ^ < |^). 



3.2.1.4 Simulation Results 

To obtain quantitative evaluation results, we consider several typical configurations 
of near-space vehicle-borne BiSAR systems. Table 3.1 gives the simulation para- 
meters with a near-space vehicle-borne receiver. The spaceborne Envisat SAR [30] 
and airborne radar are assumed as the transmitters in Case A (for transmitter A 
and receiver B) and Case B (for transmitter B and receiver B), respectively. We 
note that an imaging coverage size of dozens of square kilometers (the size is from 
2 x 3 km to 6 x 6 km in the simulation examples) can be obtained. Figure 3.6 
gives the example range resolution results. It is noted that the range resolution 
has geometry-variant characteristics, which depend on not only the slant range 
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Fig. 3.6 Range resolution of the example near-space vehicle-borne BiSAR configurations. 
a Case A Transmitter A and Receiver A. b Case B Transmitter B and Receiver B 
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Fig. 3.7 Azimuth resolution of the example near-space vehicle-borne BiSAR configurations. 
a Case A Transmitter A and Receiver A. b Case B Transmitter B and Receiver B 



but also the azimuth range. The range resolution degrades with the increase of 
azimuth range displaced from the scene center. To ensure a consistent range res- 
olution, the imaged scene coverage should be limited or a long slant range should be 
designed. 

Figure 3.7 gives the example azimuth resolution results. It is noted that the angle 
ifi') between the transmitter velocity vector and receiver velocity vector has also 
an impact on the azimuth resolution. Additionally, we note that the azimuth res- 
olution between scene center and scene edge has a small performance difference. 
This phenomenon is caused by the change in azimuth time. Further simulations 
show that, if the synthetic aperture time is determined primarily by the receiver, and 
suppose the transmitter's velocity is constant, the azimuth resolution will degrade 
with the increase in the receiver's velocity. However, if the synthetic aperture time 
is determined primarily by the transmitter, and suppose the transmitter's velocity is 
constant, the azimuth resolution will be improved with the increase in the near-space 
vehicle-borne receiver's velocity. 
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range 



Fig. 3.8 Targets with the same range delay at zero-Doppler but will have different range 
histories 



3.2.2 Azimuth-Variant Characteristics 



Unlike normal azimuth-invariant BiS AR, from Eq. 3 .27 we note that the range history 
of a point target does not depend any more on the zero-Doppler distance and the 
relative distance from the target to the transmitter, but also on the absolute distance 
to the receiver. This azimuth- variant phenomenon brings a great challenge towards 
developing efficient image formation algorithms. Consequently, the locus of the slant 
ranges at the beam center crossing times of all targets parallel to the azimuth axis 
follow an approximate hyperbola. 

We consider a typical geometry shown in Fig. 3.8 (left), where the (x, y) plane is 
locally tangent to the surface of the earth. The targets are assumed to lie on this plane, 
and the transmitter velocity vector is parallel to the y-axis. Figure 3.8 (right) shows 
the trajectories of three targets A, B, and C, as well as the hyperbolic locus. Because 
of the hyperbolic locus, after linear range cell mitigation correction the targets such 
as A, C and D are at the same range gate. This is different from the monostatic case, 
which has a linear locus instead of hyperbolic. 

Consequently, the general range Doppler image formation algorithm cannot han- 
dle this problem in a high-precision manner. Moreover, there is a range ambiguity 
that does not exist in monostatic cases. It can be observed from Fig. 3.8 that two or 
more targets (e.g., A, C and D) located at different positions can have the same range 
delay at zero-Doppler but will have different range histories (curvature). Similar phe- 
nomena have been investigated in [31], where a ground-based stationary receiver is 
assumed. For the near-space vehicle-borne BiSAR this problem becomes being more 
complex, since the spacebore or airborne transmitter follows a rectilinear trajectory, 
while the near-space vehicle-borne receiver follows also a rectilinear trajectory but 
with a different velocity. 
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The situation becomes more complicated for unfiat digital-earth model (DEM) 
topography. In monostatic or azimuth-invariant BiSAR systems, scene topography 
can be ignored in developing image formation algorithms because the measured range 
delay is related to the double target distance and the observed range curvature; hence, 
topography is only used to project the compressed image which is in slant range to 
the ground range. However, for the near-space vehicle-borne BiSAR it is mandatory 
to know both the transmitter-to-target distance and target-to-receiver distance to 
properly focus on its raw data, which clearly depend on the target height. In this 
case, the conventional imaging algorithms such as Chirp-Scaling and wavenumber- 
domain, may not be suitable to accurately focus on the collected data. 



3.2.3 Two-Dimensional Spectrum Model 

The two-dimensional spectrum model is useful for developing an effective image 
formation algorithm. We consider the demodulated near-space BiSAR signal 



S(fr, fa) = P{f, 



oo 
r) ■ / exp( - 



rt lL±A R(r) _j2 7t f aT \dr (3.34) 



cq 



where f r , f a and /o denote the range-frequency, azimuth-frequency, and carrier 
frequency, respectively. Loffeld et al. [32] divided the bistatic phase history into two 
quasi-monostatic phase histories and expended them into Taylor series around the 
individual points of stationary phase 



<D r (r) = -lit 



co 



*x2 
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(3.35) 



$ r (r) = -In 
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(3.36) 



Let 



4>;(r,*) = 0, <D;«)(T-r r *) = 



(3.37) 



There is 
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(3.38) 
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(b) 



Instantaneous transmitter squint-angle [deg.] 
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Fig. 3.9 The constraints of the Loffeld's BiSAR two-dimensional spectrum model, a In small 
squint-angle cases, b In small difference between vj and v^ 



In this way, one two-dimensional BiSAR spectrum model is constructed in [32]. 
Unfortunately, Eq. 3.38 is constrained by 



r,o! 2 « 2RJ ° 



Ivj 



r ,2 „ 2Rl rO 

Ivi 



(3.39) 



where r r o and r r rj denote the shortest transmitter-to-target distance and receiver- 
to-target distance, respectively. Figure 3.9 shows that this model can only be used 
in small squint angle or small difference between the transmitter velocity and the 
receiver velocity. 

Another model using series expansion to express azimuth time as a function 
of azimuth frequency during azimuth Fourier transform was proposed in [33], and 
a chirp scaling algorithm was developed based on this model [34]. The accuracy 
is controlled by keeping enough terms in the power series, but only the first several 
series can be used. Using the Fresnel approximation, one 2-D BiSAR spectrum model 
was derived in [35]. This model can be used only in small squint angle. Moreover, it 
has no advantages in approximate errors. In summary, most BiSAR spectrum models 
have to use some approximations by keeping enough terms in the power series, but 
cannot be used in the BiSAR with large difference between transmitter-to-target 
distance and target-to-receiver distance or large squint angle, like the near-space 
vehicle-borne azimuth-variant BiSAR. 

To develop an effective 2-D spectrum model for the near-space vehicle-borne 
BiSAR, we consider the azimuth-variant BiSAR geometry shown in Fig. 3.10 and 
the range history sum to an arbitrary point given in Eq. 3.27. Taking spaceborne 
transmitter as an example and supposing the synthetic aperture time is T s , since 
there is 



v, » v r , R t0 » v,-t, r e [-T s /2, T s/2 ] 



(3.40) 



we then have 



3.2 Near-Space Vehicles in Passive BiSAR Imaging 

^ equivalent to stationary receiver 



33 



& 



^ 



V, 



t ec^ et 




:,r^^* 




Fig. 3.10 Azimuth-variant BiSAR geometry and its equivalent model 
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Fig. 3.11 
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Range errors caused by the equivalent BiSAR geometry model, a R,q = 800 km, v, = 
bii,o = 15km, v, = lOOm/s 



R(X) = JR% + (VtT) 2 + Jr 2 q + (V,T) 2 « 7^0 + (VeqT) 2 + R r0 (3.41) 

where \ eq — \ t + \ r is the equivalent velocity between the transmitter and the rec- 
eiver. Figure 3.11 gives the range errors caused by the equivalent BiSAR model. The 
phase errors that are smaller than ?r/4 can be ignored for subsequent image formation 
processing algorithms. We thus conclude that this equivalent BiSAR model is feasible 
for near-space vehicle-borne BiSAR systems. 



3.2.4 Image Formation Processing 



The equivalent instantaneous azimuth Doppler frequency f a is represented by 



fa 



> eq 



sin (</>,) 



(3.42) 
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where <j) t is the instantaneous angle between the transmitter and the point target. As 
there is Xf a /v e q 4C 1, when 4>t is small we have 



R(f a ;R,o) 
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2 fl 
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For simplicity, we denote R(f a ', Rto) by 

R(f a ;Rto) = Rtott + C s ] + R r o 
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2v 2 
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(3.43) 



(3.44) 



(3.45) 



is the scaling factor of the chirp scaling (CS) algorithm. 

Next, a nonlinear chirp scaling (NCS) algorithm which is similar to the algorithms 
developed in [31, 34] can then be applied. Suppose the transmitted signal is 



s,(t) = exp[J7t(2f c t + k r n] 



(3.46) 



where f c and k r denote the carrier frequency and the chirp rate, respectively. Note that 
the amplitude terms are ignored. The demodulated baseband signal in the receiver is 



s r (t) = exp 



-H-^J 



Rb(r) 
co 



(3.47) 



Applying one Fourier transform (it is often implemented with fast Fourier trans- 
form (FFT)) to azimuth time r yields 



S r (t, fa) = exp - jnk 
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with 
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_L _ XR to 



(Kfa/Veq) 2 



[y/l-afa/Vc^) 2 ] 



(3.49) 



where (x p , y p , 0) is the point target's coordinate, and Rbo is the smallest bistatic 
range. 
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Chirp scaling processing with the phase term 
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co 



(3.50) 



where /?t>ref is the reference range. Next, they are transformed into two-dimensional 
frequency-domain through range FFT. After range compression with the phase term 



$ 2 = exp 



J l + C/ r 



exp 



(W,^#) (3.51) 



a range inverse FFT (IFFT) is applied. To compensate the effects of azimuth-variant 
Doppler chirp rate, phase filtering is further applied before azimuth IFFT [31] 



<t>3 = exp 
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(3.52) 



Finally, the Doppler chirp rate is corrected by 
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Applying one Fourier transform in the azimuth direction and multiplying a phase 
compensation term 
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(3.54) 



The focused BiSAR image can then be obtained by an azimuth inverse Fourier 
transform. The whole processing step is illustrated in Fig. 3.12. 

To evaluate the performance of the derived imaging algorithm, the two example 
BiSAR data from five point targets are simulated using the parameters listed in 
Table 3. 1 . A distance separation of five times the range/azimuth resolution is assumed 
in the simulation examples. Additionally, froC^ = 0) = and £,-o(t = 0) = are 
supposed. After processed by the equivalent velocity model and NCS combined 
image formation algorithm, the simulation results are given in Fig. 3.13. Note that, 
as the reference range used in the imaging algorithm is chosen to be the central target 
range, the focusing performance of the central target is better than that of the edge 
target. Even so, from the results we can conclude that the five point targets are well 
focused and the imaging performance is acceptable. 
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$i: CS phase function 



2 : range compression 




O, : phase filter 



4 : Doppler chirp rate 
correction 



<J> 5 : azimuth compression 



Fig. 3.12 Processing steps of the equivalent velocity model and NCS combined imaging algorithm 
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Fig. 3.13 Processing results with the image formation algorithm, a Case A. b Case B 



3.3 Near-Space Vehicles in Passive Environment Monitoring 



The last but not the least use of near-space vehicles in passive remote sensing is 
regional atmospheric environment monitoring. An example application is the near- 
space vehicle-borne GPS receivers for radio occultation measurements. Radio occul- 
tation measurements using GPS and a LEO spaceborne receiver can offer accurate 
profiles of atmospheric refractivity, pressure, water vapor, and temperature with a 
high vertical resolution [36, 37]. While GPS occultation data collected from space 
have the advantage of being of global coverage (one LEO spaceborne receiver pro- 
vides about 500 globally distributed occultations per day), the occultation sampling 
in any particular region is relatively sparse. In contrast, near-space vehicle-borne 
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local curvature center 



Fig. 3.14 Geometry of the near-space vehicle-borne GPS occultation 



GPS receivers can offer a dense radio occultation measurement over a specific area 
of interest for purposes of regional weather and climate observation or coupled 
ocean/atmospheric process research. 

As near-space vehicle-borne receivers operate in the atmosphere, the receivers can 
track GPS satellites as they set and rise behind the Earth's limb, therefore collecting 
GPS signals at both negative and positive elevations relative to the receiver's local 
horizon. The viewing geometry of near-space vehicle-borne GPS occultation can 
then be thought of as a hybrid between spaceborne and ground-based occultation 
geometries; hence, it combines the high vertical profiling capability of spaceborne 
occultation data (at least for height below the receiver) with the benefit of routinely 
obtaining a relatively large number of daily profiles in the regions of interest. Each 
occultation yields a profile of refractivity below the height of the receiver with a 
diffraction-limited vertical resolution. A single near-space vehicle-borne receiver 
with a full 360° field of view can observe several hundreds of occultations per 
day scattered within a radius of several hundreds of kilometers from the receiver. If 
multiple near-space vehicle-borne receivers can be formatted with an optimized flight 
plan [38], more daily occultations can be obtained over the region. This information, 
when combined with columnar water vapor distribution derived from upward looking 
ground-based receivers, and possibly moisture information from microwave remote 
sensing [39], is extremely useful for regional weather monitoring and hydrological 
and boundary layer research. 

As shown in Fig. 3.14, from the knowledge of the positions of the GPS satellite 
and near-space vehicle-borne receiver and their clocks (obtained from other GPS 
measurement collected simultaneously) the delay due to the intervening media can 
be isolated. For spaceborne GPS occultation, when the bending is determined from 
the Bouguer's formula, the index of refraction can be derived from 
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oo 

a(a') 



n{a) — exp 



it J y a /2 _ a 2 



da' 



(3.55) 



where a denotes the asymptote miss distance or impact parameter. The other para- 
meters are defined as Fig. 3.14. Since the upper limit of the integral is infinity, it is 
necessary to have measurements of a starting from outside the atmosphere (where a 
vanishes) for this integral to be performed. Differently, since the near-space vehicle- 
borne receiver lies inside the atmosphere at radius r^ , then only the measurements 
of a(fl') for a<rRii(rR) are available; therefore, in this case Eq. 3.55 should be 
modified. 

In fact, when the receiver is inside the atmosphere, the a(a') has a maximum 
exactly at the point that separates positive and negative elevation measurements (how- 
ever, spaceborne GPS occultation has no local maximum or minimum). In this case, 
the maximum value for a is a max = n(r^)r^. Then, from the Abel inversion scheme 
[40] we have 



n(a) = exp 



1 f a(a') ,' 

da 



n J V a '2 _ a 2 



(3.56) 



In this way, the numerical refractions n can be obtained. 

GPS occultation using near-space vehicle-borne receiver offers a simple method 
of obtaining vertical information about the lower atmosphere up to the receiver 
altitude. About several hundreds of measurements per day per receiver could be 
obtained in a given region, in contrast to radiosondes that are launched once or 
twice daily. Since vertical information is extremely important in characterizing the 
stability of the atmosphere and is expensive to acquire, the profiles obtainable from 
near-space vehicle-borne GPS receivers may provide quite useful applications such 
as regional weather forecasting, hydrology, surface-air exchange, and related topics 
on atmospheric research. 

However, as a new GPS occultation technique, much further work is needed. For 
example, when the near-space vehicle-bore receiver is moving inside the atmosphere, 
separating the positive and negative elevation data will be more difficult, depending 
on the motion dynamics. Additionally, to assess the merits of near-space vehicle- 
borne GPS occultation, further experiments should be performed to investigate the 
possible strategies of complementing bending measurements with refractivity and to 
understand their impact on inversion performance. 

3.4 Potential and Challenges 

Near-space passive remote sensing provides many promising potential applications 
due to the superiorities of persistent regional coverage, robust survivability, bistatic 
observation [41], and low cost. Certainly, there also are several technical challenges 
in synchronization compensation and motion compensation. 
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3.4.1 Potential Applications: Homeland Security 

One potential application is homeland security. To protect civilian population, mass 
transit, civil aviation, and critical infrastructure from terrorist attacks, significant 
improvement in our capability to safeguard homeland security is required, so that 
we can efficiently neutralize these threats without impacting the normal day-to-day 
human and commerce activities. To reach this aim, radar has long been used in a 
variety of military and civilian applications and has been an essential component 
of the current defense systems. Presently, many countries have their own civil avia- 
tion radar networks. These radar networks are specifically designed to ensure early 
warning against potential hostile targets. However, the range and nature of potential 
threat targets are becoming more and more diverse. The tragic event of 9/1 1 is but 
one example. Others might include missile attacks and UAVs. This means that the 
source of attack can emanate in a much wider variety of new and different forms, 
and cannot be efficiently dealt with by the current radar systems. It appears that pas- 
sive near-space vehicle-borne radars using opportunistic illuminators can provide a 
potential solution to these problems. 

Rather than emitting signals, this passive remote sensing system relies on oppor- 
tunistic transmitters to detect potential targets. This allows the receiver to oper- 
ate without emitting energy. This is particularly attractive for homeland security 
applications, because it is desirable for such a sensor to serve also other purposes 
such as traffic monitoring and weather prediction without influencing normal human 
activities. 



3.4.2 Potential Applications: Persistently Disaster Monitoring 

Another potential application is persistent disaster monitoring [42, 43]. The fre- 
quency of natural disasters has shown a rapid increase in the recent years [44]. 
Examples of this trend are related to floods, earthquakes, tsunamis, hurricanes, and 
forest fires [45]. The tsunami that killed thousands of people in the coastal areas of 
India, Indonesia, Thailand, and Srilanka has brought the awareness that we cannot 
regard natural disasters any more as inevitable or unavoidable and remain helpless 
observers. Methods and strategies along with effective techniques have to be devel- 
oped to predict and tackle natural disasters. Toward this aim, immediately after that 
tsunami, natural disaster monitoring has received much recognition [46], but there 
is still a lack of feasible and practicable solutions. 

Taking tsunami detection as an example, tsunamis can be triggered by earthquakes, 
submarine landslides, volcanic eruptions, meteorite impacts, or by a combination 
of all these factors. Tsunami waves have sufficient energy to cross entire oceans, 
traveling outward rapidly with small detectable heights in the deep ocean. But, since 
a tsunami wave's shoaling is near the shore, a tsunami that is visibly unnoticeable in 
the deep ocean will become more detectable if closer to the shore [47]. It is reported 
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that the Sumatra tsunami had a recorded wave height of 60-80 cm in the deep ocean, 
but had a maximum wave height of 15 m near Banda Aceh [48]. As tsunami is a 
surface gravity wave with a wavelength much larger than the ocean depth, the wave 
develops in three stages [49]: the first stage is the formation due to the initial cause 
and propagation near the source, the second is the free propagation of the wave in 
the open ocean at large depths, and the third is the propagation of the wave in the 
region of continental shelf and shallow coastal waters. The wave behavior in each 
specific region should be understood in order to properly detect, predict, and model 
the tsunami wave's activity. 

It is has been proved that, there is a link between the tsunami wave amplitude and 
the microwave RCS [44], and significant variations (a few dB) of the RCS synchro- 
nous with the sea level anomaly can be found both at C and Ku band in the geophysical 
data record of the altimetry satellite Jason- 1 . From the microwave remote sensing 
viewpoint, the oceanographic detectable features of a tsunami wave are summarized 
in [9]: (1) tsunami wave height, (2) tsunami orbital velocities, (3) tsunami-induced 
RCS modulations. The first two parameters inherently belong to an ocean wave, and 
the third is a geophysical feature that arises from complex hydrodynamic processes. 
The behavior of tsunami-induced RCS modulations depends upon bathymetry, mete- 
orological factors, and sea-state. The capability of microwave radar systems to detect 
internal waves has been proved [50-52]. Along with tides, tsunamis are shallow water 
waves and have the potential to trigger internal waves. Tsunami-induced RCS mod- 
ulations can be used to predict future tsunamis. This involves not only detection but 
also an estimate of the tsunami magnitude. Thus, near-space passive remote sensing 
has the possibility to detect a tsunami by comparing the pre and post quake patterns. 



3.4.3 Challenges: Synchronization Compensation 

In monostatic radar remote sensing, the colocated transmitter and receiver use the 
same stable local oscillator (STALO), the phase can decorrelate only over a very 
short time (about 1 x 10~ 3 sec). By contrast, for the near-space passive remote sen- 
sing, there is no phase noise cancellation because the transmitter and receiver use 
separate local oscillators. This superimposed phase noise will corrupt the received 
signal over the whole coherent processing time. 

As only local oscillator phase noise is of interest, the passive radar can be sim- 
plified into an azimuth-only system [53]. Suppose the transmitted signal is sinusoid 
with phase argument 

<D r (f) = 2nf T t + M<p T (t). (3.57) 

The first term is caused by the carrier frequency and the second term is phase noise, 
and M is the ratio of the carrier frequency to STALO frequency. After reflected from a 
target, the received signal phase is that of the transmitted signal delayed by the round- 
trip time r. The receiver output signal phase <!>(?), results from demodulating the 
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received signal with the receiver STALO which has the same form as the transmitter 
STALO 

<S> R (t) = 2nf R t + Mip R {t) (3.58) 

can be expressed as 

<D(0 = 2n(f R - f T )t + Infix + M(<p R (t) - cp T (t - r)). (3.59) 

The first term is a frequency offset arising from non-identical STALO frequen- 
cies, which will bring a drift on the final image. Because this drift can be easily 
corrected using ground calibrator, it is ignored here. The second term forms the 
usual Doppler term which should be preserved. The last term represents the effect of 
STALO frequency instability which is of interest. It is assumed that (pr(t) and <p R (t) 
are independent random variables having identical power spectral density (PSD) 
Sip(f)- Then the PSD of phase noise in a bistatic radar system is 

V(/) = 2M 2 V/) (3-60) 

where the factor 2 arises from the addition of two uncorrected but identical PSDs. 
It has been proved that synchronization compensation is required for bistatic radar 
systems [54, 55]. 

Several potential synchronization techniques or algorithms [56-58], such as using 
ultra-high-quality oscillators [59] and an appropriate bidirectional link [60] are pro- 
posed. However, in most cases, we cannot alter the transmitter; hence, it is necessary 
to develop a practical synchronization technique without alteration to the transmit- 
ters. One potential solution is the following direct-path signal-based time and phase 
synchronization approach [61]. 

Suppose the nth transmitted pulse with carrier frequency fj„ is 

x„(t) = s(t)exp(j27Tf T nt)expU<Pd(n)) (3.61) 

where (pd(n) is the original phase, and s(t) is the radar signal in baseband. Let t dn 
denote the delay time of direct-path signal, the received direct-path signal is 

s rf«( f ) = s ( f ~ tdn)exp[j2n(f T n + fdnW - tei)]exp(j <pd(n)) (3.62) 

where f dn is the Doppler frequency for the «th transmitted pulse. Suppose the demod- 
ulating signal in receiver is 

S f (t) = exp-j(2nf Rn t) (3.63) 

Hence the received signal in baseband is 

Sdn(t) = s(t - t dn )exp(-j2jt(f Tn + f dn )t dn )sxp(j2irAf n t)exp(j(p d(n) ) (3.64) 
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with Af n — fjn — fRn , here <pd(n) is the term to be extracted to compensate the phase 
synchronization errors in the reflected signals. 

Matched filtering with the reference function yields 



y dn (t) ^Af n 



sm(jrAf„(t-t dn + ^)) 

7lAf n (t-t dn + ^-\ 

- j\2rc(fdn + fRn)tdn 



exp 



nAfi 



J7tAf n (t-t dn +^\ 



<Pd(n) 



(3.65) 



A/„ 



where k r is the chirp rate of the transmitted signal . The maxima will be at t — tdn — 
where we have 



exp 



/jrA/„(f-r rf „ + ^)lL 



**-& = l 



Hence the residual phase term in Eq. 3.65 is 



*(n) = -2n(f d „ + f Rn )tdn 



Ml 



<Pd(n) 



(3.66) 



(3.67) 



As A/„ and k r are typically on the orders of 1 kHz and 1 x 10 13 Hz/s respectively 



[25], 



TtAfnl 



has negligible effects. Equation 3.67 can be simplified into 



f(n) = -lltifdn + fR„)tdn + <Pd(n) 

Similarly, we have 

f(n + 1) = -2jT(fd( n + l) + fR( n +l))td(n + l) + <Pd(n+l) 

Let 

fd(n) — fdO + Sfdn, /«(«) = //fO + &fRn 



(3.68) 



(3.69) 



(3.70) 



where fdo and /ro are the original Doppler frequency and error-free demodulating 
frequency in receiver, respectively. Correspondingly, Sfd„ and 8/r„ are the frequency 
errors for the nth pulse. Hence we have 

Win+\) ~ <Pd(n) =[f(n + 1) - if(n)} - 2jt(f R o + fdo)(td(n+l) ~ tdn) 

- 27t(Sfdn + SfR„)(td(n + l) — tdn) (3.71) 

Generally speaking, 8f dn + 8ff>„ and £</(n+l) — tdn are typically on the orders 
of 10 Hz and 10 _9 s, respectively, then 2n(&fd n + &fRn)(td(n+\) — tdn) is found to 
be smaller than 2it x 10~ 8 radian, which has negligible effects. Furthermore, since 
td( n +\) and tdn are known, Eq. 3.71 is simplified into 
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<Pd(n+l) ~ <Pd(n) = fe(n) (3.72) 

where 

fe{n) = t^C" + 1) - f(n)} ~ 27T(f R0 + fdo)(td(n + l) ~ tdn)- (3.73) 

The (pd{n) can then be obtained by iterative computation of Eq. 3.73. Next, the 
phase synchronization errors in the reflected channel can be compensated with this 
information. 

The direct-path signal-based synchronization approach provides a simple solution 
to compensate the synchronization errors, but there is a disadvantage in that the 
receiver must fly with a sufficient altitude and position to maintain a line-of-sight 
contact with the transmitter. 

3.4.4 Challenges: Motion Compensation 

In the previous discussions, we did not consider the motion errors. However, prob- 
lems will arise due to the presence of atmospheric turbulence, which introduce plat- 
form trajectory deviations from normal position, as well as altitude (roll, pitch, and 
yaw angles) [62, 63]. Consequently, motion compensation is required. In current 
radar systems, GPS and inertial navigation systems (INS) are usually employed for 
this task. By contrast, for near-space vehicle-borne radars, motion compensation 
facilities may not be reachable because near-space vehicles have very limited load 
capability. Therefore, some efficient motion error compensation techniques should 
be developed. 

For the near-space vehicle-borne BiSAR, since its synthetic aperture time is 
short, we ignore the acceleration errors in along-track and consider only the motion 
errors in cross-track in the following discussions. As shown in Fig. 3.15, sup- 
pose the ideal transmitter and receiver instantaneous positions at azimuth time r m 
are (v t r m , y t o, h t ) and (y r r m , y r o, h r ) respectively, and their actual positions are 
(v f T m , y, + Ay,(t m ), h,) and (y r r m , v r0 Ay,-(r m ), h r ). 

Suppose the transmitter motion error in cross-track is Ar t (r m ), there are 
Ay ( (r m ) = — Ar ? (r m )cos (a to) («ro is the instantaneous incidence angle from the 
transmitter to the point target P„(x n , y n , 0)), Azt(x m ) — — Ar t (r m ) sin(a r o) ar >d 



v„ - y, = r tn cos (a, ) with r tn = ^hj + (y n - y t0 ) 2 and h, = r t „ sin(a r0 )- 
The transmitter range history can then be represented by 



RT(r m ) = V (x n - v,r m ) 2 + (y n - y,o - Ay,(r m )) 2 + (h, - Az t (r m )) 2 

= yj{x n ~ v,x m ) 2 + r 2 n + 2r, n ■ Ar t (r m ) + Ar 2 (r m ). (3.74) 

Assume that the instantaneous transmitter squint angle is 6 tm and denote x t (r m ) — 
— v t r m and tan(0, m ) = x t (r m )/r tn , we can then get 
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Fig. 3.15 Illustration of the 
near-space vehicle-borne 
BiSAR motion errors 
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Similarly, for the receiver range history we have 



RR(tm) ^ J r ™ + x r( x m) ~ Ar r {x m ) COS (6 rm ) 



(3.75) 



(3.76) 



where r m , x r (r m ), Ar r (r m ) and 6 rm are defined in a manner similar to r tn , x t (x m ), 
Ar r (r m ) and 9 tm . The bistatic range history can then be represented by 



R(r m ) 



, + xf(x m ) - Ar,(x m ) cos (0, m ) 
+ yjr^ + x?(x m ) - Ar,\x m ) cos (0 nn ). 



(3.77) 



As the first and third terms are the ideal range history for subsequent image 
formation processing, we consider only the second and fourth terms. We have 

d[Ar,(x m ) cos (6 tm ) + Ar r (r m ) cos (0 m )] 



Ar,(r m i) 
Ar r (r m i) 



dx„ 



sin(9 tm ) ■ v, ■ (T m2 - T m i) 



sm(6 rm ) ■ v r ■ (x m2 - T m i) 



(3.78) 



As there are Ar t (x m \)/ r tn <K 1 and Ar r {x m \)/r rn <K 1, when x m i — x m \ is short, 
Eq. 3.78 will be equal to zero. That is to say, the motion errors in a short interval 
can be seen as constant. This phenomenon validates again the equivalent velocity 
and NCS combined imaging algorithm. Otherwise, efficient motion compensation 
techniques should be applied. A potential solution is the subaperture-based motion 
compensation algorithms. The details can be found in [64, 65]. 
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3.4.5 Challenges: Antenna Directing Synchronization 
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Another technological challenge is antenna directing synchronization, which requires 
the transmit antenna and receive antenna to simultaneously illuminate the same 
region on the ground. Although the feasibility of the BiSAR concept has already 
been demonstrated by experimental investigations [66], BiSAR antenna directing 
synchronization is still a technological challenge. 

The concept of pulse chasing was proposed as a means for bistatic radar antenna 
directing synchronization [67, 68]. In this approach, the single receive beam rapidly 
scans the area in two dimensions covered by the transmit beam, essentially chasing 
the pulse as it propagates from the transmitter. Because the single beam chases one 
pulse at a time, this imposes a limit on the maximum allowable pulse repetition fre- 
quency (PRF). Moreover, antenna directing synchronization is even more demanding 
for BiSAR systems, especially for azimuth-variant BiSAR systems, since the geo- 
metric and radiometric characteristics of azimuth-variant BiSAR are strictly dictated 
by illuminator configuration and operation. The requirement of antenna directing 
synchronization for interferometric Cartwheel SAR was simply analyzed in [69]. 
Several altitude and antenna directing synchronization strategies for satellite forma- 
tion configuration were proposed in [70]. The comparison of altitude and antenna 
pointing design strategies of noncooperative spaceborne BiSAR were investigated 
in [71]. 

Without loss of generality, we consider a rather general BiSAR configuration, in 
which the transmitter and receiver are mounted on different platforms. As a typical 
example, suppose the normalized transmit antenna gain is 



Gt(0t, <Pt) = exp 



I"/ e T \ 2 / <f>T \ 2 1] 

L\#r3dB/ K&TMB/ J J 



^r3dB/ \0T3dB, 

where y is the Gauss parameter which is often assumed to be 1.3836, 9j and 4>t are 
the antenna beamwidth in range and azimuth respectively. Correspondingly, 9jidB 
and 4>tmb are their 3dB beamwidth. If there are range antenna directing errors A9j , 
we then have 

AGt(0t,<Pt) 4y6 T A0 T 



Gt(0t,4>t) 

Similarly, for the receive antenna we have 
AGr(0r,<Pr) 



Then, there is [72] 
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Fig. 3.16 Zoomed impact 
of antenna directing 
synchronization errors 
in range 
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To evaluate the impact of antenna directing synchronization errors on SAR imag- 
ing performance, linear and quadratic errors are usually assumed in existent papers. 
In fact, antenna directing error usually is oscillatory for practical systems. Hence, 
we use an oscillatory model 



A0 T {t) = A6 R (t) = A r cos (cM + <Po) 



(3.83) 



with A r , <D r and cpo are the amplitude, frequency and initialization angle in range, 
respectively. Take an azimuth-invariant X-band BiSAR, in which the transmitter and 
receiver are moving in parallel tracks with a constant and identical velocities v s — 
150m/s using the following parameters: 6r(t) — 8°, 0R(t) — 5°, A r — 1, &>,• = 
1 .98 and c/>o = as an example, the corresponding impacts are illustrated in Fig. 3.16. 
It is seen that antenna directing synchronization errors manifest themselves as a 
deterioration of the impulse response function. They may defocus BiSAR image and 
introduce a significant increase of the sidelobes. 

Similarly, suppose the azimuth antenna directing synchronization errors are rep- 
resented by 



Afi T = &Pr = A a cos (co a t + fa) 



(3.84) 



where A a , a> a and fio are the amplitude, frequency and initialization angle in azimuth, 
respectively. The corresponding transmit/receive azimuth antenna figure can be rep- 
resented by 



w(t) = Wait - A a sin(oV + A))] 



(3.85) 



with w a {t) — sinc 2 (^^(f — t c )), where L a , R c and t c are the antenna length, 
nearest range from the platform to the scene and its corresponding time, respectively. 
From the SAR processing procedure, we know that the azimuth signal after range 
and azimuth compressions is 
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T s /2 

S out (t) = expO-2,/^ - jntf) . / „,<T)capC*rA*T)dt (3.86) 

-T s /2 

Since 

r » A a sm(co a t + Po) (3-87) 

According to the principle of Taylor series expansion, we have 

wit) = w a (t) - A a sin(<w fl f + /So) ■ u)„(0 (3.88) 

Substitute Eq. 3.88 to Eq. 3.86 we get [73] 

T s/2 

Soutit) =expij2jtf c t - J7tk r t 2 ) ■ / w a ir) expi2jik r tr)dr 



-T s /2 

-^iJlnf c t-jnkS).A a J sin( Wfl r + ^)exp(2^rr)dr. 

-T,/2 

(3.89) 
We can see that, similar to range antenna directing synchronization errors, azimuth 

antenna directing synchronization errors also introduce paired echoes. 

Thus, a high-precision antenna directing synchronization approach should be 

developed for near-space vehicle-borne azimuth-variant BiSAR systems. 

3.5 Conclusion 

In this chapter, we considered the near-space vehicle-borne passive radar in sur- 
veillance and reconnaissance. It has the advantages of entire planet coverage and 
simple transmitter-receiver synchronization, but high-resolution performance can- 
not be obtained. Thus, we further considered near-space vehicle-borne BiSAR high- 
resolution imaging. A comprehensive description of the system imaging performance 
was provided, along with the signal processing challenges. Also, we considered near- 
space vehicles in a passive environment monitoring. Potential applications in home- 
land security and disaster monitoring and technical challenges in synchronization 
compensation and motion compensation are investigated. 
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Chapter 4 

Near-Space Vehicles in High-Resolution 

Wide-Swath Remote Sensing 



Abstract Spaceborne SAR has an imaging capability of wide- swath (the width of 
the ground area covered by the radar beam) with a limited azimuth resolution. By 
contrast, airborne SAR has an imaging capability of high azimuth resolution, but 
limited swath coverage. There is, therefore, a desire to increase swath coverage and 
azimuth resolution simultaneously. As near-space vehicles operate at altitudes higher 
than that of airplanes but lower than satellites with a high flying speed, compared 
to spaceborne and airborne SARs, simultaneous relative high-resolution and wide- 
swath (HRWS) remote sensing is possible for near-space vehicle-borne SAR. In this 
chapter, we explained how near-space vehicles could be exploited for future HRWS 
remote sensing applications. 

Keywords Near-space • High-resolution wide-swath • Multiple apertures • Wave- 
form diversity • Multiple-input and multiple-output (MIMO) • Space-time coding 

4.1 Restrictions on Achievable Resolution and Swath 

An efficient remote sensing technique should provide high resolution imagery over 
a wide area of surveillance, although there is a contradiction between azimuth reso- 
lution and swath width [1, 2]. A good azimuth resolution requires a short antenna to 
illuminate a long synthetic aperture which results in a wide Doppler bandwidth. This 
calls for a high pulse repeated frequency (PRF) to sample the Doppler spectrum. 
Thus, azimuth resolution and ambiguity suppression impose a lower bound on the 
PRF, and the higher it is, the better the achievable azimuth resolution becomes and 
the smaller the ambiguous signals will be. But, a high PRF means a smaller swath 
width. The relationship between the maximum imaging swath on ground, W s , and 
the required PRF can be expressed as 

en 

W s < (4.1) 

2 ■ sin(»7) ■ PRF 

where cq and r] denote the speed of light and the incidence angle, respectively. 
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Substituting the expression of SAR azimuth resolution p a and rearranging the 
terms in Eq. 4. 1, we then have 

^ < C ° (4.2) 

p a 2 ■ v s ■ sin(^) 

where v s is the SAR platform velocity. Generally, cq/v s is nearly constant at 20,000 
for LEO spaceborne SARs and typically in the range of 300,000-750,000 for air- 
borne SARs. Near-space vehicles can fly at a speed ranging from stationary to 
1,500 m/s, the corresponding cq/v s will be greater than 100,000. Thus, compared 
to spaceborne and airborne SARs, the near-space vehicle-borne SAR provides a 
more flexible choice between azimuth resolution and swath width for HRWS remote 
sensing. 

Equation 4.2 can also be reformed into the basic minimum antenna area 
constraint [3] 

4v s XR c ■ tan(^) 

^antenna == "a ' L a > (4.3) 

CO 

where H a is the antenna width, L a is the antenna length, X is the radar wavelength, 
and R c is the slant range from the radar to the mid-swath. This requirement arises 
because the illuminated area of the ground must be restricted so that the radar does 
not receive ambiguous returns in range dimension or/and azimuth dimension. In 
this respect, a high operating PRF is desirable to suppress azimuth ambiguity. But 
the magnitude of the operating PRF is limited by range ambiguity requirement. 
Therefore, in conventional SARs the unambiguous swath width and the achiev- 
able azimuth resolution impose contradicting requirements on system design, and 
consequently allows only for a concession between azimuth resolution and swath 
width. These considerations can be combined into one inequality which expresses 
the range of slant range values that are appropriate to each of the unambiguous swath 
intervals 

(£ +I >)!^(w-'')? 

where n a is an integer that effectively labels each swath interval, and T p is the pulse 
duration. 

Although near-space vehicle-borne SAR provides a more flexible choice between 
azimuth resolution and swath width, the illuminated swath must be restricted so that 
the received radar echoes are unambiguous in range dimension or/and azimuth dimen- 
sion. As an example, assuming an X-band near-space vehicle-borne SAR system with 
the following typical parameters: X — 0.03 m, r] — 60°, v s = 1,000 m/s, L a — 
0.4 m, H a — 0.22 m, and the flying altitude is Ho — 60km (according to [4], 
Ho — 60 km is the best altitude for high-speed near-space vehicles), then the maxi- 
mal swath width is found to be 32.7 km. 



4.1 Restrictions on Achievable Resolution and Swath 53 

The attainment of a wider swath will become increasingly difficult if a higher 
azimuth resolution is required, due to the requirements of increased PRE The PRF 
limit for a certain swath W s and a given incidence angle r\ is [5] 

2v ? en 

— -<PRF<0.S - (4.5) 

L a 2VK. r sin(??) 



4.2 State-of-the-Art: HRWS Remote Sensing 

Numerous innovative concepts are presented to enable HRWS remote sensing 
whereof the most promising ones employ a single transmit antenna in combina- 
tion with multiple channels on receive [5-1 1]. This means that these systems have a 
receive antenna which is split into multiple subapertures with independent receiver 
chains that are interpreted as individual channels. A similar technique is the distribu- 
tion of the receiver apertures on multiple platforms that are grouped into a formation 
of space or airborne sensors leading to a multistatic SAR [12-15]. It is common 
to all of these methods that the backscattered signals are received simultaneously 
by multiple apertures, which are mutually displaced in azimuth and/or elevation 
dimension. The different receive aperture positions introduce spatial diversity in the 
received echoes. The basic idea is, hence, to use the multiple receivers to gather 
additional information and to benefit from this information to overcome the above 
restrictions of conventional SARs [16-23]. In general, the multiple receivers can be 
either arranged in flight direction ("along-track"), perpendicular to it ("cross-track"), 
or in both dimensions. They are illustrated respectively in Fig. 4.1. 



4.2.1 Multiple Apertures in Elevation 

The multiple apertures in elevation proposed by Griffith and Mancini [16] consists 
of an array antenna split in elevation. The overall antenna dimension is smaller 
than that implied by the minimum antenna area constraint yielding a broad beam in 
elevation dimension that covers a wide swath but at the same time gives rise to range 
ambiguous echoes. The range ambiguities are suppressed by adaptively steering nulls 
in the antenna pattern in elevation to the directions of the ambiguous returns. In this 
way, a widened swath can be obtained. However, in a monostatic SAR the swath will 
be no longer continuous because blind ranges are introduced when the receiver is 
switched off during transmission. This restriction can be overcome by using bistatic 
SAR which allows for simultaneous transmission and reception [24-26]. 

In 2001, Suess et al. [5] proposed an innovative processing approach in elevation 
dimension for HRWS remote sensing. The SAR system is built on an array antenna 
consisting of multiple elements in elevation. The processing concept combines the 
echoes from the subapertures in elevation in a way to form a narrow beam in elevation 
which "scans" the ground in real-time in order to follow the echo of the transmitted 
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Fig. 4.1 Typical configurations of the multiple apertures or receivers 



signal on ground. This enables the suppression of range-ambiguous returns and 
ensures a high antenna gain. However, this technique requires knowledge of the 
observed terrain topography, or else a mispointing of the narrow elevation beam may 
occur, resulting in severe gain loss [27]. 



4.2.2 Multiple Channels in Azimuth 



In 1992, Currie and Brown [2] proposed the displaced phase center antenna (DPCA) 
in azimuth technique, which is based on dividing the receive antenna in along-track 
direction into multiple subapertures, each receiving, down-converting, and digitiz- 
ing the radar echoes. Hence, for every transmitted pulse the system receives multiple 
pulses in the along-track direction. This means that additional samples can be gath- 
ered, thus increasing the effective sampling rate on receive. Consequently, either the 
resolution can be improved while the swath width remains constant, or the swath 
can be widened without increasing azimuth ambiguities or impairing the resolution. 
That is to say, the system benefits from the whole antenna length regarding azimuth 
ambiguity suppression, while azimuth resolution is determined by the dimension of a 
single subaperture, thus decoupling the restrictions on HRWS remote sensing. How- 
ever, this approach imposes a stringent timing requirement on the system regarding 
the relation between sensor velocity, PRF, and antenna length. 
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Regarding multiple channels in azimuth, the approach presented in [28, 29] intro- 
duces a phase correction that is applied to the raw data to resample the signal in 
azimuth. This method is based on an analysis of the multichannel signal's phase 
compared to the phase of a monostatic and uniformly sampled signal. This yields a 
Doppler frequency-dependent phase difference between the multichannel signals and 
the monostatic signal. Hence, by applying an appropriate phase correction to the data 
of each individual channel, the phase of the multichannel signal is adjusted in such 
a way that the resulting phase corresponds to the monostatic and uniformly sampled 
signal. In 2003, Krieger et al. [30] proposed an algorithm based on a generaliza- 
tion of the sampling theorem that allows for unambiguous recovery of the azimuth 
spectrum from multiple aliased subaperture signals. This method is elaborated in 
several follow-on papers [31-36] and extended to the burst operation mode [37]. 

4.2.3 Multiple Apertures in Two Dimensions 

In 1999, Callaghan and Longstaff [17] proposed the quad array approach which is 
based on an antenna split into two rows and two columns yielding a four-element 
array. This approach can be understood as a combination of the two approaches 
described previously. Thus, this approach combines the advantages of gathering 
additional samples in azimuth to suppress azimuth ambiguities and simultaneously 
enabling an enlarged swath for a fixed PRE However, the proposed system in eleva- 
tion will also result in blind ranges in the imaged swath. Additionally, the stringent 
timing constraint to ensure a uniform spatial distribution of the gathered sampled in 
azimuth dimension must also be satisfied. 

Besides, Claassen and Eckerman [38] proposed an alternative concept based on 
steering multiple beams to different azimuth directions and assigning each of the 
corresponding footprints to a different slant range. In other words, the footprint of 
each beam steered to a different squint angle corresponds to a different subswath 
of the overall imaged region. As the antenna elevation dimension is larger than the 
imaged swath, range ambiguities can be well suppressed, thus allowing for a PRF 
high enough to suppress the azimuth ambiguities. However, there will be coarsened 
resolution and impaired performance arising from the needed high squint angles [22]. 

4.2.4 Distributed SAR Constellations 

The concept of sparse array SAR constellations was proposed in [12, 13, 23]. This 
approach considers an arbitrarily distributed spaceborne array of radar satellites that 
do not necessarily build a formation aligned in the along-track dimension. The satel- 
lites are then regarded as a sparse antenna array that collect also angle-of-arrival 
information. This additional information enables wide-swath imaging with enhanced 
azimuth resolution. However, the processing of the irregularly spaced and sparsely 
distributed samples of the array requires a space-time minimum mean-squared error 
estimator. An optimum way of processing in the space-time domain is derived 
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Fig. 4.2 Geometry mode of the near-space vehicle-borne SAR for HRWS remote sensing 

in [13]. Further, Aguttes [14] proposed the concept of SAR train consisting of a 
multi-satellite constellation distributed in along-track direction and uses a spread 
spectrum waveform for transmission. 



4.3 Near-Space Vehicle-Borne SAR HRWS Remote Sensing 

To reduce the required PRF for near-space vehicle-borne SAR HRWS remote sensing, 
we consider first the multiple azimuth beam technique. The incentive is that it allows 
an effective decrease in the operating PRF of a SAR system, while ensuring that 
the bandwidth of target signals is still adequately sampled in azimuth. Two different 
approaches (see Fig. 4.2) to implementing this technique are described below. 

4.3.1 Single-Phase Center Multibeam SAR Imaging 



The operation mode of the single phase center multibeam (SPCM) SAR system is to 
transmit pulses by a single broad azimuth beam and receive the returns by a number 
of narrow contiguous azimuth beams that span the mainlobe width of the transmit 
beam, as shown in Fig. 4.3. A distinct channel is associated with each of the receive 
beams, and hence, the data are split according to the azimuth angular position or, 
equivalently, instantaneous Doppler frequency center in the azimuth direction. As 
a result, given the knowledge of the relative squint angles of each beam (hence the 
Doppler center frequency for each beam) and assuming suitable isolation between the 
beams, each channel can be sampled at a Nyquist rate appropriate to the bandwidth 
covered by each narrow beam, instead of that covered by the full beamwidth. 

This arrangement enables correct sampling of the azimuth spectrum with a PRF 
fitting the total antenna azimuth length, which is N times smaller than the PRF 
necessary for the antenna azimuth length, L a . The area of each beam antenna is then 
restricted by 
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Fig. 4.3 Geometry mode of 
the near-space vehicle-borne 
SAR for HRWS remote 
sensing. (Reproduced by 
permission of © 2009 
Elsevier Masson SAS) 
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Clearly, the minimum area of each beam antenna is A^-times smaller than the respec- 
tive area of a monostatic SAR. Correspondingly, the relationship expressed in 
Eq. 4.2 is changed into 



W K 



Nc 
2 ■ v s ■ sin(?7) 



(4.7) 



From Eq. 4.7 we note that the relation not only depends on the platform velocity v s 
and the incidence angle r\ but also on the number of antenna beams. 

Thereafter, the DPCA technique discussed in [39, 40] can be used to gain addi- 
tional samples along the synthetic aperture which enables an efficient suppression of 
the azimuth ambiguities, i.e., the multiple beams in azimuth allow for the division of a 
broad Doppler spectrum into multiple narrow-band subspectra with different Doppler 
centroids. A coherent combination of the subspectra will then yield a broad Doppler 
spectrum for high azimuth resolution. Thus this technique is especially attractive for 
high-resolution SAR that uses a long antenna for unambiguous wide-swath coverage. 

The data in each channel will be aliased to the relative zero-frequency. This 
problem can be resolved by using the prior knowledge of the relevant Doppler center 
frequencies to regain the full azimuth Doppler bandwidth. Consider Fig. 4.3 [41] 
which shows a near-space vehicle-borne SAR with three beams in azimuth. For a 
given coordinate of the SAR, the imaged area is simultaneously illuminated by the 
three beams. We may view this SAR as a conventional SAR (the central beam), 
operating with a PRF which is one third of that required to adequately sample its 
beamwidth, together with two additional beams on either side of the central one. As 
described previously, the basic idea is that the additional samples obtained by the 
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outer beams will fill in the gaps in a target's azimuth phase history which occur due 
to operating at a low PRE 

For the left beam, its instantaneous slant range in terms of the slow time variable 
r is given by [42] 

V T COS 

Rl + v 2 r 2 - 2R v s x sin0 s ^ R Q - v s r sin<9, + -*-— (4.8) 

2ao 

where 7?o lS the nearest slant range of the target, and 6 S is the squint angle. The 
Doppler centroid frequency and azimuth bandwidth can be derived from Eq. 4.8 as 

2 2v s cos 5 
fdi = -v s sin6 s , B d i = — (4.9) 

where L as denotes the length of each beam antenna. Similarly, for the central beam 
and right beam, we have 

(4.10) 
(4.11) 



Generally, we have the following approximations because the 6 S is small. 

Bdl ^ Bdc ^ Bd\ — Bds, fell - fdc — fdc - fdr — Bds (4.12) 

The ambiguous Doppler spectrum can be recovered unambiguously by applying a 
system of reconstruction filters. Some reconstruction algorithms have been proposed 
by other authors, e.g., [34, 31]. A block diagram for the reconstruction from the 
three-channel signals is shown in Fig. 4.4. This algorithm is based on considering 
the multichannel near-space vehicle-borne SAR data acquisition as a linear system 
with multiple channels, each described by a linear filter. From the sampling theorem, 
we know that the sampled signal spectrum X s (/) is the sum of the unsampled signal 
spectrum, Xo(f), it repeats every f s Hz with f s the sampling frequency. 

There is 

00 

X s (f) = fsY J X (f-nf s ) (4.13) 



If fs > 2Bds, the replicated spectra will not be overlapped, and the original spec- 
trum can be regenerated by chopping X s (f) off above f s /2. Thus, Xo(f) can be 
reproduced from X s (f) through an ideal low-pass filter that has a cutoff frequency 
of fs/2. The subsequent data processing involves interpolating the data from each 
channel. We then have 



ys(n) 



U). n=0,±N ,±ZN ,... (4U) 

0, otherwise 
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Fig. 4.4 Multichannel 
signal reconstruction 
algorithm in case of three 
channels 
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where 2Vo is the interpolation scaling factor. Accordingly, the spectrum is represented 
by 



W)= Z ys(n)e- j2nfn = X s (N f) 



(4.15) 



Then, the linear filter is derived as 



H k (f) 



No, 
0. 



/ 



fck 



<^,*€(U,r) 



(4.16) 



otherwise 



where A: e (/, c, r) is shown in Fig. 4.4, and f c k is the corresponding Doppler 
centroid. 

Finally, the filtered signals can be combined coherently, as shown in Fig. 4.5. 
That is to say, in a manner similar to adaptive antenna beamforming technique, the 
filtered signals can be combined coherently. In this way, the capability of ambiguity 
suppression allowing for improved resolution and an enlarged swath can be achieved. 
Note that, for optimum performance the relationship between sensor velocity and 
along-track offsets of the three subchannels must result in equally spaced effective 
phase centers , so that a uniform sampling of the received signal can be obtained [31]. 



4.3.2 Multiple Phase Center Multibeam SAR Imaging 



The multiple phase center multibeam (MPCM) SAR system also synthesizes multiple 
receive beams in the azimuth direction; however, the operating mode of this system 
is quite different from that of the previous one. In this case, the system transmits 
a single broad beam and receives the radar returns in multiple beams which are 
displaced in the along-track direction. The motivation is that multiple independent 
sets of target returns are obtained for each transmitted pulse if the distance between 
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Fig. 4.5 Azimuthal spectra synthesis for multichannel subsampling in case of three channels. 
(Reproduced by permission of © 2009 Elsevier Masson SAS) 



phase centers is suitably set. This method basically implies that we may broaden the 
azimuth beam from the diffraction-limited width, giving rise to improved resolution, 
without having to increase the system operating PRE 



4.3.2.1 Imaging Scheme and Signal Model 

Consider Fig. 4.3, for the central beam there is 



Re(x) = jRn+ v 2 Ax - T ) 2 - 2R Q v s {r - r ) sin( 



Rq - v s sin# s (r - r ) 



v 2 cos 2 S 



2Rn 



-(.T-ToY 



(4.17) 



where To is the azimuth reference time. The phase history can then be represented 

by 



*c(T) 



AttR c (x) 

X 



(4.18) 



Suppose the spatial separation between two phase centers is d a , the phase histories 
of the left beam and right beam can be represented, respectively, by 
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n th pulse position 
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Fig. 4.6 Equivalent phase centers for the case of three receive beams 
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(4.20) 



From the DPCA principle, we known that subsequent data processing simply 
involves interleaving the data from each stream properly in the azimuth direction 
so that the proper relative phasing is maintained. Azimuth compression can then 
be processed as normal. This technique may be extended to more than three receive 
beams, the main effect being that the constant phase difference increases in magnitude 
for the additional beams. Note that the application of this technique assumes that 
cubic (and high order) terms may be neglected in the azimuth phase history of targets. 
A general criterion for this condition to hold is that this term causes less than a jt/2 
excursion in phase over the aperture synthesis time. 

To construct the synthetic aperture, the system must operate with a PRF which 
leads to a properly sampled synthetic aperture appropriate to the beamwidth of the 
system. In fact, the operating PRF of the system is always equal to the Nyquist rate 
for the diffraction-limited beamwidth of the antenna, regardless of the number of 
receive beams and the width of the beams. This is illustrated in Fig. 4.6 for the case 
of three receive beams. 



4.3.2.2 Non-uniform Displaced Phase Center Sampling 

Taking into account that the overall antenna length in azimuth L a is made up of N 
apertures, each length is d a , the optimum PRF can be derived from 

2v 5 2v s 
PRFuni = T7^ = -j 1 (4.21) 

Nd a L a 

This imposes a stringent requirement on the system as it states that to ensure equal 
spacing between all samples in azimuth the PRF has to be chosen such that the SAR 
platform moves just one half of its antenna length L a between subsequent radar 
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pulses. This optimum PRF yields a data array equivalent to that of a single-aperture 
system operating with N ■ PRF. In reverse, any deviation from the relation will 
result in a non-equally sampled data array along the synthetic aperture that is no 
longer equivalent to a monostatic signal and cannot be processed by conventional 
monostatic algorithms without performance degradation. 

To analyze the impact of nonuniform displaced phase center sampling, we con- 
sider the received radar returns 



st(t, r) » ai[ho(t) ®, h it i(t, t)], i = 1,2, 



N 



(4.22) 



where a, is the RCS parameter, t is the range fast time, r is the azimuth slow time, 
and®, is a convolution operator on the variable t. ho(t) and/zi ,(?, r) denote, respec- 
tively, the range reference function and azimuth reference function 



ho(t) — w r (t) ■ exp I — jnk r t 2 \ 



(4.23) 
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In 



x w fl (r) -8 



R c (r) + R c (r + i^\\\ 
Rc(r) + R c (r + ij) 



CO 



(4.24) 



where w r {t) and w a (r) denote the antenna pattern in range dimension and azimuth 
dimension, respectively. As there is 

J? c (t) + Rc (r + i — J « 2R C (r + ij-\ (4.25) 



we then have 



si(t, r) « a 



ho(t) ® t h 



(— M 



(4.26) 



with 
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{'' r + l t) = Wa ( r + l t) • exp {- J T Rc ( T 



2v ? 



co 
Equivalently, the nonuniform PRF can be considered as azimuth time drift. 

da 1 



(4.27) 
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Fig. 4.7 Reconstruction filtering for multichannel subsampling in case of three channels 
After matched filtering and range mitigation correction, we get 
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where k is an integer, fy is the Doppler frequency centroid, and k a is the Doppler chirp 
rate. It is noted that the signals are periodic nonuniform with the period of l/PRF. 
This information is particularly important for developing nonuniform reconstruction 
algorithms. 



4.3.2.3 Azimuth Signal Reconstruction Processing 

The ambiguous Doppler spectrum of a nonuniformly sampled SAR signal can be 
recovered unambiguously by applying a system of reconstruction filters. The algo- 
rithm illustrated in Fig. 4.7 is based on considering the data acquisition in a DPCA 
SAR as a linear system with multiple receiver channels, each described by a lin- 
ear filter. The reconstruction consists essentially of multiple linear filters which are 
individually applied to the subsampled signals of the receiver channels and then 
combined coherently. 

Taking three azimuth beams as an example, the corresponding reconstruction 
filters are 
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where T = d a /(2v s ), K = exp(j2jtT ■ PRF) - 1, and G = expC/"4jrT ■ PRF) - 1. 
The reconstructed signals can be expressed as 

Sa \wrf) = Zl Ismf) ® PI te) + Z2 (sprF) P2 \3PRf) 

+ z 3 ( )®p 3 ( ) (4-33) 
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(4.34) 

(4.35) 
(4.36) 



where pi (■), P2O) an d P3(') are me time-domains representing Pi (/), P%(f) and 
Pl(f), respectively. 



4.3.3 Ambiguity -to-Signal Ratio Analysis 



For a given range and azimuth antenna pattern, the PRF must be selected such that 
the total ambiguity noise contribution is adequately small relative to the signal, i.e., 
from —18 to — 20 dB. A low PRF will increase the azimuth ambiguity level due to 
increased aliasing of the azimuth spectra. On the other hand, a high PRF value will 
reduce the interpulse period and result in overlap between the received pulses in time. 
To resolve these problems, the transmit interference restriction on the PRF must be 
satisfied with 
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n> n' + 1 

< PRF < -^ — (4.37) 



2 -K min T ^R 



P CO 

where R m m is the nearest slant range, R max is the farthest slant range, n' is a given 
integer, and T p is the pulse duration. 

Similarly, the nadir interference restriction on the PRF can be represented by 

< PRF < ^ =j- (4.38) 
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where m' is a given integer and h s is the platform altitude. Alternatively, given a PRF 
or range of PRFs, the antenna dimensions and/or weighting (to lower the sidelobe 
energy) must be adequately small such that the ambiguity-to-signal ratio (ASR) 
specification is satisfied. 

The azimuth ambiguities arise from finite sampling of the Doppler spectrum at 
intervals of the PRF. Since azimuth spectrum repeats at the PRF intervals, the signal 
components outside this frequency interval will fold back into the main part of the 
spectrum, and the desired signal band will be contaminated by the ambiguous signals 
from adjacent spectra. This can be evaluated by the azimuth ambiguity to signal ratio 
(AASR) defined as [43] 

Z~=-oo J_°o. 5 % G 2 (f + m ■ PRF) df 

AASR « ^ ^ (4.39) 

f-™ d Bc GHf)df 

where B ( ] and G(f) denote the SAR correlator's azimuth processing bandwidth and 
azimuth antenna pattern, respectively. 

Considering the three beams illustrated in Fig. 4.3, we have 



( nL as 



Jl l^ ns cos(( ■ 6*) \ 
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(4.40) 

Note that i e (—1,0, 1) is determined from the positions of the three beams, i.e., 
i = — 1 is for the left beam, i = is for the central beam and i = 1 is for the right 
beam. Because the 3dB beamwidth is approximately determined by [43] 

0*^f, B s =k a ^- (4.41) 

with k a a given constant. Equation 4.40 can then be further simplified into 
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Fig. 4.8 The AASR of an example near-space vehicle-borne SAR as a function of PRF. (Repro- 
duced by permission of © 2009 Elsevier Masson SAS) 



Then, from Eq. 4.39 we can get 
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(4.43) 
An example near-space wide-swath SAR is considered with the following para- 
meters: h s = 60 km, v s = 1, 000 m/s, k — 0.03 m, L as — 1.2 m, the calculated 
AASR results are illustrated in Fig. 4.8. Note that the AASR is typically specified 
to be on the order of — 20 dB; however, even at this value ambiguous signals may 
be observed in images that have very bright targets adjacent to dark targets, because 
SAR imagery can have an extremely wide dynamic range due to the correlation 
compression gain. As such, a lower AASR, e.g., — 30 dB, is desirable. 

Similarly, range ambiguities may result from the preceding and succeeding pulse 
echoes arriving at the antenna simultaneously with the desired return. This type of 
ambiguity is relatively not significant for near-space SAR. The range ambiguity of 
the near-space vehicle-borne HRWS SAR can be analyzed similar to the general 
SAR, which has been fully investigated, e.g., [3]. 
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Thus the set of PRFs is established by the acceptable maximum range and AASR 
requirement, as well as the transmit and nadir interference. Note that, there may be 
no acceptable PRFs at some incidence angles that meet the minimum requirements. 
The designers then have the option to relax the performance specifications for the 
imaging area or exclude these modes from the operations plan. 



4.3.4 Conceptual System Design 

An example system is designed to manifest the imaging performance. It operates in 
X-band with a center frequency of 10 GHz. From the radar equation, we obtain [24] 

NESZ = s — '-' / (4.44) 

Pa Vg NAl s p r 

where/? r is the average slant range which is assumed constant, K (K « 1.38 x 10 -23 ) 
is the Boltzmann constant, T sys is the system noise temperature, L f is the loss factor, 
F„ is the receiver noise figure, P avg is the average transmit power, N is also the 
number of subapertures, A as is the sub-aperture antenna area, and p r is the range 
resolution cell for one look. 

To calculate the system performance we assume a range resolution of p r —0.2 m, 
an overall loss factor of L f — 3 dB, and a receiver noise figure of F„ — 3 dB. It is 
further assumed that the signal bandwidth is adjusted for varying incidence angles 
such that the ground-range resolution is constant across the whole swath. An example 
system design is provided in Table 4.1. We note that, for the incidence angle given 
in Table 4.1, the swath width is about 8km and the NESZ is approximately — 48 dB. 
These results show that a satisfactory performance can be achieved, however, with 
only a small number of subapertures which have relatively small antenna area. The 
conclusion is that unambiguous range and swath width can be obtained using a 
near-space vehicle-borne SAR with multiapertures in azimuth. 



4.4 Near-Space HRWS Remote Sensing via Multiple Apertures 

The DPCA technique allows either the azimuth resolution to be improved while 
the swath width remains constant, or the PRF to be reduced without increasing 
the azimuth ambiguities. However, it imposes a stringent timing requirement on the 
system regarding the relationship among platform velocity, PRF, and antenna length; 
otherwise, the azimuth signal will be nonuniformly sampled. As digital beam forming 
(DBF) on receive is a promising candidate for HRWS remote sensing [44-47], in 
this section we describe one reflector antenna-based DBF solution to near-space 
vehicle-borne SAR HRWS remote sensing. 
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Table 4.1 Performance 
parameters of an example 
near-space vehicle-borne 
SAR 



Parameters 



Variables 



Values 



Mean transmit power 


Pavg 


1W 


Number of sub-aperture 


N 


3 


Sub-aperture antenna length 


i- 1 as 


0.40 m 


Sub-aperture antenna width 


H a 


0.10m 


Near-space vehicle velocity 


Vs 


500 m/s 


Incidence angle 


n 


20° 


Swath width 


w s 


6.80 km 


Radiometric resolution 


NESZ 


-48.58 dB 


Incidence angle 


n 


25° 


Swath width 


w s 


7.30 km 


Radiometric resolution 


NESZ 


-48.11 dB 


Incidence angle 


n 


35° 


Swath width 


w s 


8.94 km 


Radiometric resolution 


NESZ 


-46.73 dB 




in azimuth 






swath Doppler spectra 

Fig. 4.9 Reflector antenna based system architecture with DBF in elevation and azimuth 



4.4.1 System Architecture and Imaging Scheme 

In SAR remote sensing missions the antenna is a key element in the total system 
performance. Designing one customized antenna with features such as wide swath 
and high resolution is a technical challenge in current SAR systems. As SAR antenna 
acts like a two-dimensional spatial filter, requirements imposed on both azimuth 
pattern and elevation pattern are important. To obtain HRWS imaging for the near- 
space vehicle-borne SAR, a reflector antenna with multiple feed elements and digital 
beamforming technique are employed. As shown in Fig. 4.9, the reflector antenna 
based near-space vehicle-borne SAR involves digital beamforming in elevation and 
digital beamforming in azimuth which are discussed, respectively, in the following 
two subsections. 



4.4.1.1 Digital Beamforming in Elevation 

To implement wide-swath imaging, the reflector antenna based digital beamforming 
in elevation is employed. The reflector antenna consists of a parabolic reflector and 
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a feed array of transmit/receive elements. To illuminate a given angular segment in 
elevation, only the corresponding feed elements are activated. On transmit activating 
all elements generate a wide beam illuminating the complete swath. On receive 
the reflected signals illuminate the entire reflector antenna, but they are focused 
on individual feed elements. The receive beam scans the complete swath within 
the time period of one pulse repetition interval l/PRF, whereas each element is 
only active during a subinterval of this time period. When a high PRF is employed, 
multiple portions of the swath will be illuminated instantaneously but each actives 
a different subset of the feed elements because of different angle-of-arrival [48, 49]. 
Equivalently, a wide swath can be obtained. 

For a strictly spherical Earth model, the angle-of-arrival of a point target echo is 
associated by 

,, . f (h s + R e ) 2 + R 2 e + R 2 (td) \ ,,,,, 

"^ = arCC ° S ( 2(h s + Re) R( td) ) (4 - 45) 

where R e is the Earth radius and R(t c i) with t c \ the two-way time delay is the slant- 
range distance to one given target. Thus, there is a one-to-one relation between 
the required beam steering angle r\ s {fd) an d the time variable fy. This information 
provides a potential solution to wide-swath SAR imaging. It is well known that, to 
avoid range ambiguities resulting from the preceding and succeeding pulse echoes 
arriving at the antenna simultaneously, the slant range R should be [50] 



2 



I— '— + T P + AT tr ) < R < — ( — T„) (4.46) 

\PRF ' ) 2 \PRF p ) 



where ; is an integer, T p is the pulse duration, and AT tr is the switching time between 
the transmission and reception of a pulse. 

Consider a SAR which operates with a PRF appropriate to the desired azimuth 
resolution requirement, but with a wider (M times) swath than that implied by 
Eq. 4.46, we then have 



co / i 
2 \PRF 



\ en (i +M \ 

T p + AT lr y R< -(—-T p } (4.47) 



As the idea of scan-on-receive (SCORE) based digital beamforming technique is to 
shape a time varying elevation beam in reception such that it follows the echo of the 
pulse on the ground [5], we divide the whole swath into M sub swaths. 
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As a compromise between computation complexity and imaging performance, 
here M is determined by 

2tfmax ■ PRF 2R miR ■ PRF 
M = j- max .J _ L mm j ^ ^ 

co co 
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where [■] and |_-J denote respectively the maximum and minimum integers, R maK 
and R m i n denote the maximum and minimum slant ranges within the imaged swath. 
To avoid range ambiguity, for the conventional SAR system the PRF should be 
satisfied with 

n' n' + 1 

< PRF < -= (4.50) 



*MTiin T 1 T -"-max i r 

cq 1 P Tr P c ~*~ Tr P 



where n' is a given integer, and x tp is the receiver protecting window extension about 
the pulse duration, T p . However, by applying the SCORE operation mode, for the 
reflector antenna-based SAR system the PRF can only be satisfied with the relation 

ri ri + 1 

< PRF < -= (4.51) 



Rm T ~ R, 



T„ - T, 



7)1 + 1 
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where R m and R m +\ are the slant ranges to two adjacent subswaths, respectively. 
Comparing Eqs. 4.50 and 4.5 1 , we can see that, for the same operating PRF a higher 
Doppler bandwidth can be sampled allowing for an improved azimuth resolution 
while keeping the range ambiguities constant. Alternately, the PRF can be reduced 
without an increase of azimuth ambiguities and degradation of the azimuth resolution 
while increasing the unambiguous imaging swath width. 

Suppose the looking-down angles of the first formed subaperture to each subswath 
areai(r), <X2(r), . . . , a,M(r), the relative phase delay from the first subswath to each 
formed subaperture can then be represented by [51] 

2nd r sm(a\(r)) 2it(M - \)d r sin(ai(r)) 

x '■■■' X { } 

where d r denotes the elevation distance between two subapertures. Similarly, for the 
second subswath we have 

2-ndr sin(o!2(r)) 2it(M - \)d r sin(a 2 (r)) 

U, , . . . , ^T-.JJJ 

A A 

Similar relations can be obtained for the remaining subswaths. They can be formed 
as a matrix expressed in Eq. 4.54, 



A r {r) 



exp 



/■ antrrinfri) ) exp ^M|*llj exp(; 2jrrf ' s i n(aM) ) 



exp ^j 2nW-l)d r sHm) j exp | j 2»(M- l H i i,fe) j £xp ^ 2jr i M = M^n( ml \ 

(4.54) 

where 

/ en / i + m — l\\ 
a m = a\J-h+ — \y l<m<M (4.55) 
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with 



a(x) 



2x(h s + R e ) 



(4.56) 



where / represents the sampling time (fast time). As an analog-to-digital conver- 
ter (ADC) is placed after each T/R-module in the feed array, a posteriori, digital 
beamforming on receive can then be formed in the direction of a wanted subswath. 
Equivalently, a large swath can be synthesized. 



4.4.1.2 Digital Beamforming in Azimuth 

To further alleviate the requirements of HRWS imaging imposed on the minimum 
antenna area, digital beamforming in azimuth is further employed. In planar antenna- 
based digital beamforming in azimuth systems, all subapertures cover the same angu- 
lar segment, thus " seeing" the identical Doppler spectra. Consequently, for the same 
PRF a higher (e.g., N times) Doppler bandwidth can be sampled allowing for an 
improved azimuth resolution while keeping the range ambiguities constant. Alter- 
natively, the PRF can be reduced by 1/7V while increasing the unambiguous range 
by a factor of TV thus extending the equivalent imaging swath. However, from the 
sampling theorem, we know that the sampled signal spectrum is the sum of the 
unsampled signal spectrum. 

As the Doppler spectra is undersampled, the Doppler spectrum can only be recov- 
ered unambiguously through the combination of the total spatial samples. That is 
to say, subsequent azimuth processing must combine the total N receive channels, 
each subsampled with PRF and aliased in frequency domain, to a single channel of 
N ■ PRF that is free of aliasing. Moreover, the relationship between platform veloc- 
ity and along-track offsets of the azimuth subchannels must result in equally spaced 
effective phase centers, otherwise, there will be nonuniform spatial sampling which 
makes the subsequent Doppler spectrum synthesis a challenge. 

In contrast to conventional planar antenna-based systems, for the reflector antenna- 
based digital beamforming in azimuth systems, there are single transmit feed and 
multiple receive feeds that are displaced in the along-track direction. Each azimuth 
element illuminates at a different angle and covers a distinct angular segment. Thus, 
each element samples a narrow Doppler spectrum. To exploit the large antenna array 
for signal transmission, the multidimensional waveform encoding investigated in 
[52] is modified and applied in azimuth direction. As shown in Fig. 4.10, in azimuth 
a series of subpulses instead of a wide duration pulse are transmitted, each subpulse is 
separately transmitted with a different short time delay by different transmit beams. 

Note that the PRF must be high enough such that the spatial sampling of each 
beam or channel is adequate. If the Doppler spectra of the elements are contiguous, 
they can jointly yield a higher azimuth resolution, D/(2N), where D is the reflector 
antenna diameter. 
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Fig. 4.10 Modified 
multidimensional waveform 
encoding in azimuth 



azimuth 




4.4.2 Imaging Performance Analysis 

As discussed previously, the antenna area will affect the azimuth resolution and 
the available swath width. Additionally, antenna beam shape, especially its sidelobe 
characteristics, is also key to the imaging performance. Ambiguity noise is thus 
an important consideration. The constraints expressed in Eqs. 4.46 and 4.47 are 
approximate and the required signal to ambiguity noise ratio may not be met. Thus, 
it is necessary to analyze the RASR and AASR performances. 



4.4.2.1 RASR Performance Analysis 



At a given time t within the data record window, range ambiguous signals arrive from 
the ranges of 



en ( m + / — 1 \ 
Ri = — ( — — + t\,m = l,2,...,M;j = ±l,±2,.., 



±N h (4.57) 



where j, the pulse number (J — for the desired pulse), is positive for the preceding 
pulses and negative for the succeeding ones, j — N/, is the number of pulses to 
the horizon. We can see that this expression is different from the expression for 
single-aperture SAR systems, which is 



Ri 



2 \PRF 



±l,±2,...,±N h . 



(4.58) 



The RASR is determined by summing all signal components within the data record 
window arising from the preceding and succeeding pulses, and taking the ratio of 
this sum to the integrated signal return from the desired pulse 



4.4 Near-Space HRWS Remote Sensing via Multiple Apertures 

U 



RASR 



. -IRmm Ar 1 R 3 jSin(c/j) 






Gpg Q 



73 



(4.59) 



min Rq sin(o:o) 



where yj m is the beamforming gain of the mth subswath at the range of Rj, Gj is 
the cross-track antenna pattern at the yth time interval of the data recording win- 
dow at a given cuy , Oj is the corresponding normalized backscatter coefficient. The 
Go, <7o, Rq and «o are the corresponding parameters of the desired unambiguous 
return. 



AASR k (PRF) 
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n=—oo 
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+z 
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(4.60) 



4.4.2.2 AASR Performance Analysis 



The desired azimuth signal will also be contaminated by the ambiguous signals 
coming from adjacent spectra. For the reflector antenna-based SAR system, the 
AASR for the £fh beam can be derived as Eq. 4.60 [43], where (k, j) e [1, 2, 3, . .. , ] 
and 



Gkij(f) 



j. ( L„ ■ cos(ri k /j) ■ (f - fdck/jT 
\ 2vs 



(4.61) 



is the azimuth antenna pattern, L a is the azimuth antenna length, r\ k n and fdc,k/j are 
the squint angle and Doppler frequency centroid for the k ory'th channel, respectively. 
As the antenna steering angle is computed using the correspondence between the 
signal delay and angle of arrival, any steering error will result in a SCORE loss. 
An antenna pattern loss factor (Fk, j) is thus employed in Eq. 4.60, which can be 
calculated as [48] 



I\./ = 



r(j-k)Td+T ps /2 
J(j-k)z d -T ps /2 



jT ps 



/2 



W2 



\Gr (r - g) \ 2 d 



k*J 



(4.62) 
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where x c ] is the time delay between two adjoining subpulses, T ps is the subpulse 
duration, and G r (r) is the receive antenna's elevation pattern. 



4.4.2.3 NESZ Performance Analysis 

The target reflected power available at the near-space vehicle-borne receiver antenna 
is determined by 

_ p, ■ Gtim) op x 2 ■ G r (m) ,. 63) 

' ' 4 • Tt ■ B*(rn) ' 4 • jt ■ Agfa) ' 4 ■ * 

where P, is the transmit peak power, G ? (??,•) and G r (r)j) with ??, incidence angle are 
the transmit and receive antenna gain respectively, and Ro(r]i) is the slant range. As 
the total data samples are processed coherently to produce a single imaging resolution 
cell, the receiver thermal noise samples can be taken as independent from sample to 
sample within each pulse, and from pulse to pulse. After coherent range and azimuth 
compression, the final image SNR can be represented by 

SMW = P t -G t (r, l) .G r (r ll) .tf.c .T p .PRF.o* 

g 256 ■ ;r 3 ■ Rfai) ■ v s ■ wKxu) ■ K ■ T sys ■ B„ ■ F n ■ L f 

The NESZ can then be calculated by 

256 ■ 7T 3 ■ Rn(ni) ■ v s ■ sin(?i,) • Kq ■ T svs ■ B n ■ F n ■ L s 
NESZ = ^-^ — , — — -. (4.65) 

p, -Gtim) -G r (m) -^- co- T P -PRF 



4.4.3 Conceptual Examples and Simulation Results 

To evaluate the quantitative performance of near-space vehicle-borne SAR with 
reflector antenna for HRWS remote sensing, we consider an example system. The 
SAR operates in X-band with a carrier frequency of 10 GHz. Table 4.2 gives the 
corresponding system parameters. Note that, a far-field, flat-earth, free-space, and 
single polarization model is also assumed. It is also assumed that the SAR moves at a 
constant velocity and operates in stripmap mode. On transmit activating all elements 
gives a wide low gain beam illuminating the complete swath. On receive the energy 
returned from a narrow portion of the ground illuminates the entire reflector, but is 
focused on individual elements of the feed aperture because the SCORE technique 
is employed in the elevation. 

It is worthwhile to compare the RASR performance of the reflector antenna and 
SCORE technique combined SAR to conventional single beam SAR. Using the sys- 
tem parameters listed in Table 4.2, Fig. 4.11 shows the comparative RASR perfor- 
mance as a function of slant range. Note that equal parameters of PRF (= 4,000 Hz) 
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Table 4.2 Near-space 
vehicle-borne reflector 
antenna SAR system 
parameters 



Parameters 


Values 


Units 


Carrier frequency 


10 


GHz 


Platform velocity 


500 


m/s 


Platform altitude 


20 


km 


Reflector antenna length 


1.5 


m 


Reflector antenna width 


0.8 


m 


Earth radius 


6370 


km 


Minimum slant range 


40 


km 


Maximum slant range 


100 


km 


Transmit peak power 


1000 


W 


Pulse duration time 


10 


|xs 


Pulse bandwidth 


100 


MHz 


Range sampling frequency 


150 


MHz 


Number of azimuth channel 


3 


- 


Number of subswath 


4 


- 



and subswath width are assumed in the simulations. For efficient SAR imaging, 
the RASR should be smaller than — 20 dB. For example, the calculated RASR are 
—35 dB for the reflector antenna and SCORE technique combined case and — 12 dB 
for general single-beam case, respectively. These results clearly show that a sig- 
nificant RASR performance improvement is obtained for the reflector antenna and 
SCORE technique combined approach. This means that a wider swath is possible for 
the near-space vehicle-borne SAR with reflector antenna, however, without decreas- 
ing the operating PRF which means that a higher azimuth resolution can be obtained. 

It is also worthwhile to compare the AASR performance of the reflector antenna 
SAR with the conventional single azimuth beam SAR. Consider the system para- 
meters listed in Table 4.2, Fig. 4.12 gives the corresponding comparative AASR 
performance as a function of PRF. Note that r^ = — 10 dB is assumed in the simula- 
tion examples. In SAR remote sensing applications, AASR is typically specified to be 
on the order of —20 dB, but a lower AASR is desired. It can be noted from Fig. 4.12 
that the AASR is typically below —20 dB with a low operating PRF requirement. 
This means that a wider swath can be obtained. 

Another imaging performance is theNESZ, which is given in Fig. 4.13. Note that 
L s — 3 dB and F„ — 3 dB are assumed in the simulation. The notches in the curve 
shape are caused by the switching of the active elements and the corresponding 
pattern switching. We can see that, for the system parameters given in Table 4.2 
the swath width ranges from — 54 dB to — 44 dB. When compared to the reflector 
antenna-based spaceborne SAR investigated in [48], here the NESZ performance is 
improved by 20 dB. The reason is that the near-space vehicles are 10-20 times closer 
to the targets than LEO satellites. 

Unlike conventional planar antenna SARs in which complex multiplication and 
summation are required to form a time varying beam, digital threshold detectors are 
employed in the reflector antenna SAR because the reflectivity variation in different 
scenes causes a variation of the average power level at the receiver [48] . The threshold 
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Fig. 4.11 Comparative RASR results between SCORE on receive and conventional single beam 
as a function of slant range 
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Fig. 4.12 Comparative AASR results conventional single azimuth beam SAR as a function of PRF 



should ensure that at each time instance only the information relevant channels are 
summed up. However, in this case the azimuth Doppler spectrum will be undulated 
by the azimuth antenna pattern, as shown in Fig. 4.14. Consequently, the imaging 
performance will be degraded. This problem can be resolved by inverse filtering the 
azimuth multichannel signals. 
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Fig. 4.13 NESZ results as a 
function of slant range 
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Fig. 4.14 Impact of azimuth 
antenna pattern on reflector 
antenna SAR Doppler 
spectra 
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4.5 Near-Space HRWS Remote Sensing via Waveform Diversity 



Recently, researches on multiple-input and multiple-output (MIMO) SAR have been 
drawing more and more attention [53-56]. MIMO is a technique used previously 
in communications to increase data throughout and link range without additional 
bandwidth or transmit power. Given that MIMO SAR is in its infancy, there is no clear 
definition of what it is. It is generally assumed that independent signals are transmitted 
through different antennas, and that these signals, after propagating through the 
environment, are received by multiple antennas. Generally speaking, MIMO SAR 
has two advantages when compared to the traditional SARs. The first advantage lies 
in spatial diversity gain and flexible spatial transmit beampattern design. The second 
advantage is resolution improvement. These are the reasons why MIMO SAR system 
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Fig. 4.15 Different 
combinations between two 
chirp signals: a k n = k n , 

U K r , = k n 
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and the associated exploitation methodology can address the shortcomings of current 
SAR systems for wide-swath remote sensing applications. 



4.5.1 Waveform Diversity Design 



Unlike conventional phased array SAR systems, in M1MO SAR systems each antenna 
transmits a unique waveform, orthogonal to the waveforms transmitted by other 
antennas. As linearly frequency modulated (LFM) waveform has been widely utilized 
in SAR remote sensing due to its good properties such as high range resolution, 
constant modulus, good Doppler tolerance, and implementation simplicity, from a 
practical point of view we think that M1MO SAR should use LFM-based waveforms, 
so as to reduce the tough requirement of high transmit power. An adaptive LFM 
waveform diversity was proposed in [57], only the simple up- and down-chirp signals 
are allowed. It can be extended into orthogonal frequency diversion multiplexing 
(OFDM) chirp waveform. 

From a practical point of view, suppose the chirp signals have equal frequency 
bandwidth and inverse or equal chirp rate (i.e., \k r i\ = \k r 2\), there are four dif- 
ferent combinations between two chirp signals, as shown in Fig. 4.16. For MIMO 
SAR range compression, the auto-correlation between identical chirp signals is used 
to evaluate the range resolution, and the cross-correlation between different chirp 
signals is crucial for suppressing the range ambiguity. 

As an example, Fig. 4.16 shows the correlation characteristics of the different 
chirp combinations shown in Fig. 4.15. Note that the parameters of bandwidth B — 
10 MHz, starting frequency f s % — Hz and pulse duration T p — 5 



lfr 6 s are 
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Fig. 4.16 Correlation characteristics of the different chirp combinations: a and b have equal chirp 
rates, c and d have inverse chirp rate 



assumed in the simulations. For MIMO SAR imaging, a quantity of importance is 
the relative level between the correlation of identical chirp signals and the correlation 
of different chirp signals. Hence, the chirp waveforms with high cross-correlational 
suppression are desired. The correlation results for different starting frequencies and 
equal chirp rate are shown in Figs. 4. 16a and 4. 16b. It is obvious that the performance 
of cross-correlation suppression improves with the increase the separation between 
two starting frequencies (i.e. f s j — f s j). The correlation results for different starting 
frequencies and inverse chirp rate are shown in Figs. 4.16c and 4.16d. It is also 
obvious that the maximum and occupied time of the correlation will decrease with 
the increase of the separation between two starting frequencies. 

It is deduced from Fig. 4.16 that using chirp signals with equal chirp rate and 
I fsi — fsj I > 25 offers a satisfactory suppression of the cross-correlation compo- 
nents. However, for a specific requirement of frequency bandwidth, this also means a 
wider total transmit/receive bandwidth for the MIMO SAR radio frequency hardware 
systems. Fortunately, the chirp signals with inverse chirp rate and adjacent starting 
frequency (see Fig. 4.17d f s i — f s j — 2B) can provide a good cross-correlational 
suppression. Therefore, it is possible to use an arbitrary number of chirp signals 
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Fig. 4.17 Example OFDM chirp waveforms. The number of the chirp signals is for illustration 
only 



which occupy adjacent starting frequency (i.e., f s , — f s j — 2B) for inverse chirp 
rate or non-overlapping frequency bands. Figure 4.17 illustrates two example OFDM 
chirp waveforms. 

They can be expressed as 



r~if/ i\ 

s^t) = J]«(0 = X rect y ' exp j2jt V Skt + 2 krkt ) 

fc=l k=[ L P-l I \ ' 



(4.66) 



where f Sk and k n denote the starting frequency and chirp rate for the k chirp signal, 
respectively. We use the ambiguity function of the OFDM chirp waveforms defined 
in [58] 



K K ~ 



s k (t)s*(t + t d )eJ 2n ^dt 



(4.67) 



where /^ denotes the Doppler frequency. This equation can be calculated by numer- 
ical integral. 

Consider the OFDM chirp waveforms shown in Fig. 4. 17 and suppose the follow- 



ing parameters: B — 10 MHz, T p — 10/Lts, f S] 



70 MHz, f^ 



70 MHz, f n 



50 MHz, f S2 
10 MHz, f S4 = 30 MHz, f S5 = 50 MHz, f S6 = 10 MHz, f Sl 
30 MHz, Fig. 4.18 shows the ambiguity function characteristics. As the value of 
X(0, 0) represents the matched filtering output without any mismatch, the sharper 
the function | x (t , f c i)\, the better the range and azimuth (Doppler) resolution that can 
be obtained for the radar system. The results show that the OFDM chirp waveform 
has a satisfactory ambiguity function performance in range resolution and Doppler 
frequency resolution. This is particularly valuable for the subsequent MIMO SAR 
image formation processing. 
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Fig. 4.18 Ambiguity functions of the OFDM chirp waveforms shown in Fig. 4.17 



4.5.2 MIMO SAR-Based Wide-Swath Remote Sensing 

The multiple antennas in the MIMO SAR systems can be arranged in along-track 
(azimuth, flight direction), cross-track (elevation, perpendicular to along-track), or 
in both dimensions. In the following discussions the multiple antennas are placed 
only in the cross-track direction. If each antenna transmits a sufficiently wideband 
signal, each transmit antenna and receive antenna pair can then act as an individual 
SAR, and the MIMO SAR can then serve as multiple SARs operating independently. 
The cumulative wide bandwidth can offer an extremely high range resolution. This 
architecture has the potential to increase azimuth resolution while enabling high gain 
on target. 



4.5.2.1 System Scheme 

As shown in Fig. 4.19a, the basic idea is to form multiple transmit and receive 
beams which are steered toward different subswaths. This is different from a phased 
array radar because orthogonal transmit signals are employed in the MIMO SAR 
system. Multiple transmit and receive beams can then be formed by grouping the 
array elements into multiple groups, each forming a transmit or receive subaperture. 
These subapertures can be disjointed or overlapped in space, but disjoint subapertures 
are assumed. Suppose the transmit array has M collocated elements, whereas the 
receive array has Af collocated elements. A total of N t < M transmit beams can 
then be formed. The waveforms used in any two subapertures are orthogonal, but 
the same waveform is used in each subaperture. As shown in Fig. 4.19b (shows 
only two beams, it is for illustration only), N t transmit subapertures are used to 
form N t different directional beams, each modulated by a different and orthogonal 
function Sk(t),k — 1, 2, . . . , N t . As the signal transmitted by the Mi subaperture 
steering toward an angle 6q is given by bf(6o)sk(t) with £>&(•) the steering vector, 
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Fig. 4.19 Geometry of near-space vehicle-borne MIMO SAR with multi-antenna in elevation 



the equivalent baseband signal seen at a location with angle 6/ can be expressed as 
bl (0i)bt(do)sk(t), where (-) T denotes a transpose operator. Correspondingly, the 
overall signal is a superposition of the transmissions from all N t subapertures 



N, 



S N 



(t) = ^bl(ei)b* k (e )sk(t) 



(4.68) 



t=i 



In this way, multiple virtual transmit beams can be formed. 
As the transmitted waveforms are orthogonal, the receivers can form N r separate 
receive patterns in elevation, each with nulls in the directions of N r — 1 ambiguities, 
and a beam in the direction of the appropriate echo. Like the multiaperture SAR dis- 
cussed previously, this means that we can divide the whole swath into N r subswaths 



2 



PRF 



AT tr \ 



c (i +k 

< Ri < — | 

2 \PRF 



4 



<k < N r 



(4.69) 



For M transmit array elements and N receive array elements, there will beMxiV 
different echoes, each is weighted by a different amplitude and phase. Each receive 
channel signal is separately demodulated, digitized, and stored. Next, a posteriori, 
digital beamforming on receive can then be carried out by a joint spatiotemporal 
processing. In this way, multiple pairs of virtual transmit and receive beams can be 
formed simultaneously in the direction of a wanted subswath. 



4.5.2.2 Matched Filtering and Multibeam Forming 

Multibeam forming in elevation must be carried out for subsequent image formation 
processing. The beamforming filter coefficients should be optimally designed, so 
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Fig. 4.20 The scheme of MIMO SAR multi-beam forming in elevation 



that the desired beam direction to the intended subswath can be obtained. To reach 
this aim, the multi-beam forming scheme shown in Fig. 4.20 is used. 

Suppose the transmit array is a uniform linear array, the transmit signals expressed 
in Eq. 4.68 can be reformed into 



s(?) = vf((p t ) [si(t), s 2 (t),..., s Nt (t)Y 



(4.70) 



where v t (</>,) = [1, e ~ WW ***<!*)/*■, . . . , e -WN t ~l)lWsin<t> t )/Uf with ^ the inci . 
dence angle is the transmit array response vector, N t is the number of the transmit 
subapertures, [s\(t), S2(t), ■ ■ ■ , SN,(t)] is the signals transmitted by the transmitter. 
The receive signals can be written as a vector x(f) = [jci(0>*2(0> • ■ ■ > x N r (t)] T (N r 
is the number of the receive subapertures) 



x(0 = v r (<Pr)v( (<Pt)s(t - r d ) 



(4.71) 



where \ r ((p r ) = [1, e-^K^WA, . . . , e -;27rW-i)[(dsin0 r )A]-|T and tj represent 
respectively the receive array response vector and the time it takes the signal to travel 
the transmitter-target-receiver distance. 

The matched filtering of each receive channel can be represented by 

z nr n, it) = x nr (0 ® s* t (t) = F~ [ {F{x nr {t)} ■ F{s* t (0}} (4.72) 

where n t e [1,2,..., Nt], n r e [1,2,..., N r ], the <S>, F~ x and F denote respec- 
tively the convolution operation, inverse Fourier transform operation, and Fourier 
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transform operation. The term z n ,.n, if) is the matched filtering output at the receive 
subaperture n r and matched to the orthogonal waveform transmitted by the transmit 
subaperture n,. The total matched filtering results can then be represented in vector 
form as 

sm[nB{t - t d )] ., . T/ , . 

z r (f) = A a — — v r Qp r )v t J (</>,) 

jtB(t -t d ) 



diagle j2,Tfli '-' d \e j - 



2jzf 2 (t-t d ) J27Tf Nt (t-t d )\ 

I (4.73) 



where A a denotes the amplitude term, fi(i = 1, 2, . . . , Ni) is the starting frequency 
of each subchirp signal. 

Next, after beamforming on transmit processing, the equivalent beam can be 
represented by 

sin[jrZ?(f — td)] ,">_/ t, , i 

Zin(0 = z r (0 ■ MM = A a ^—^ ^y r (ct>r)e' 2nUn(, - td) 

JtB(t -t d ) 

sin[N t Bjt(t - td)] 

x — — V — (4.74) 

sin[5jr(f - td)] 

where \t(<pt) is the transmit steering vector and / m ; n is the smallest frequency center 
among the N t orthogonal subchirp signals. Correspondingly, for the nth receive 
subaperture we can get 

Zout,n(0 = v£ n (</v) ■ z in (f) 

smlirBit — td)] sm[N t B7T(t — td)] ;-)_/■ ,, ,\ 

= N r A a — — — — ■ e ]27lf ™" (, - ,d) (4.75) 

nB(t-td) sm[Bir(t - td)] 

where v r , n (0r) is the steering vector for the nth receive subaperture. Then there is 

sm[jrB(t - td)] sin[N t B7T(t - td)] 



Zout,n(0 



N r A a - 



nB{t — t d ) sin[B7r(f - td)] 



(4.76) 



Therefore, if the matched filtering outputs are fed into multiple digital filters 
with their coefficients designed for different beam-pointing directions, multiple pairs 
of virtual transmit and receive beams in predetermined directions can be formed 
simultaneously. Moreover, high range resolution is also obtained. We conclude also 
that a smaller ( 1 /L r for one subaperture and \/N r for the whole aperture) RCS target 
can be detected for the MIMO SAR. Correspondingly, for a given requirement of 
S AR image SNR, a lower peak transmit power or average transmit power is required. 



4.5.2.3 Numerical Simulation Results 

According to multi-transmission and multi-reception in the elevation scheme 
described previously, eight subchirp signals are transmitted from the eight transmit 
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Fig. 4.21 The returned 
signal spectrum received by 
a subaperture 
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subarrays simultaneously and received by four subapertures. The reflected signals 
will be received by all the receive subantennas. 

For simplicity, suppose there is a point target at the swath center and with an 
azimuth distance of m, Fig. 4.21 shows the signal spectrum received by a single 
subantenna. It can be noted that, after beamforming on transmit, an echo with wide 
frequency width can then be received by each receive subantenna. As the transmit sub- 
chirp signals are orthogonal, their matched filtering can be implemented separately 
with the corresponding reference functions. In this way, the final range resolution can 
be approximately improved by a factor of N t , i.e., the number of the subchirp sig- 
nals transmitted simultaneously, as shown in Fig. 4.17. As an example, suppose the 
inter-element spacing is half wavelength, Fig. 4.22 shows the magnitude of the direc- 
tional gain for an example subaperture with L t — 8 elements. Next, one or multiple 
beams can be obtained by cooperative digital beamforming on receive. Thereafter, 
the subswath data can be coherently combined into the raw data for unambiguous 
wide-swath imaging. 



4.5.3 Space-Time Coding MIMO SARfor 
High-Resolution Imaging 



In this section, we further proposed a space-time coding MIMO SAR for high- 
resolution remote sensing. As shown in Fig. 4.23, the approach employs MIMO 
configuration in elevation direction and a space-time coding scheme in azimuth 
direction, along with OFDM waveform diversity and DPCA techniques. The basic 
idea is to divide the total transmit antenna elements into multiple groups, each 
forming a subbeam, thus offering several benefits, including improved ambigu- 
ity suppression for wide-swath imaging, improved SNR, and flexible operational 
configuration. 
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pulse compression results before beamforming on transmit subaperture 
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Fig. 4.22 The final matched filtering results for one subaperture: the above one is before beam- 
forming on transmit subaperture, the below one is after beamforming on transmit subaperture 



Fig. 4.23 Slidelook 
space-time coding M1MO 
SAR (The labeling a...f 
and 1 ... 6 are just two 
different along-track position 
of the same 6-element array 
detailed in Sect. 4.5.3.3) 




4.5.3.1 Space-Time Coding Transmission in the Azimuth Direction 



The simplest Alamouti code is suitable for the space-time coding MIMO-OFDM 
SAR system with two transmit antennas [59], which is based on the following 2x2 
orthogonal matrix 
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S(Jl,S2) = 



Si 
*2 



(4.77) 



where (•)* denotes the complex conjugate operator, s\ and s 2 denote the two trans- 
mitted signals. The entry at the ith row and yth column represents the signal to be 
transmitted from the fth antenna at the y'th time interval. That is to say, two dif- 
ferent signals si and s 2 are transmitted simultaneously from the antennas 1 and 2 
respectively during the first signal period, following which the signals s% and — s* 
are transmitted from antennas 1 and 2 respectively. We assume that the target's RCS 
remains constant over consecutive signal periods. The channel response for the «th 
receive antenna can then be represented by 



H= [h\ <n , h 2 , n ] 



(4.78) 



where h\ n and h 2n denote the channel response of the two transmitted signals, 
respectively. 

The signals received by the nth receive antenna over consecutive signal 
periods are 



r\n — S\® hl >n + S 2 <S> «2,n + «l,n 
r%n — Sj <E> «l,n — S* ® «2,n + «2,n 



(4.79) 



where n\ „ and n% n denote individual additive white Gaussian noise. This signal 
model is formulated in time domain. Transforming them into frequency domain 
yields 



(4.80) 



'Rln' 




'Si s 2 ~ 




'#!,„" 




'N hn ~ 


_Rln_ 




. °2 °1 _ 




. ff 2,« . 


+ 


_N 2 ,n_ 



where R ln (R 2n ), Si (S 2 ), S*_ (S x *), H Un (H 2 , n ) and N Un (N 2 , n ) denote the 
Fourier transforming representation of n« (r 2 n), s i fe), -?| ( s *)> h\ >n («2,n) 
and ni„ (n 2t „), respectively. Next, the received radar echoes should be separated 
by decoding processing. As the transmitted signal matrix is known for both the 
transmitter and the receiver, the decoding matrix can be easily constructed as 



D = 



S 2 

-Si 



(4.81) 



Since Si and S 2 are orthogonal, the decoded signals can be represented by 



'Rln' 




~s* S 2 ' 




'Rln' 




"(ISil 2 


_Rln' _ 




_S*-Si_ 




Rln _ 




.(ISil 2 



\S 2 \ 2 )H hn 
\S 2 \ 2 )H 2 , n 



Ni, 
Ni, 



- N 2 , n 
+ N 2 „ 



(4.82) 

It is noted from Eq. 4.82 that, after Alamouti decoding, the two received signals 
have the same transmit information (|.Si| 2 + IS2I 2 ) and dependent channel response. 
It can also be concluded that the Alamouti scheme extracts a diversity order of 2 if 
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two orthogonal signals are transmitted with the same power means |,Si| 2 = IS2I 2 . 
After matched filtering, the reconstructed R[„' and Rin achieve a gain of 12 dB 
while the noise level w 6dB. Thus, the equivalent diversity gain is 6dB. 



4.5.3.2 MIMO Configuration in the Elevation Direction 

In the elevation direction, the total antenna array elements are divided into multiple 
subarrays. These subarrays can be disjoint or overlapping in space, but here disjoint 
subarrays are assumed. Orthogonal waveforms are transmitted from different subar- 
rays to illuminate a wide swath. Each receive subarray can receive all the reflected 
signals. If N t transmit subarrays and N r receive subarrays are employed, there will 
be N t x N r different returns for the receiver. For the mth (m — 1 , 2, . . . , N t ) transmit 
subarray and nth (n — 1, 2 N r ) receive subarray, the radar response or scatter- 
ing function can be approximated as a realization of one random process denoted by 
a n>m (m — 1,2,..., N t ; n = 1,2,..., N r ). The random process is such that over 
the time period of the transmit signal duration, the measured amplitude should be 
the same, i.e. |a„, m | is infact a constant (for one set of simultaneous measurements 
employing different codes). 

The received signal at the nth receive subarray due to the mth transmit waveform 
can be represented by 

r n ,m (t) = s m (t, n)a n , m + "«,m 0) (4.83) 

The term s m (t, n) denotes the signals transmitted by the mth subarray and received 
by the nth subarray. The term n nm (?) is an additive noise process being independent 
of the response function a n , m . As there are N t transmit signals, the received signals 
at the nth receive subarray are the linear combinations of all such signals 

N t N, N t 

r n = ^ r n>m (t) = ^ s m (t, n)a n , m + ^ n n>m {t) = s(n)a„ + n„ (4.84) 

m=\ m=l m — \ 

The terms s(n), a„ and n„ are the matrices whose columns are the s m (t, n), a, um 
and n, um it), respectively. As there are a total of N r receive subarrays, the data from 
each subarray can be composed into a single vector 

R = rf,rJ,...,r£ r 7 '=S-A + N (4.85) 

where T denotes the transpose, S = diag[s(l), s(2), . . . , s(A^)], A = [ai, &2, . . . , 
aN r ] T and N = [ni, n2, . . . , njy,.]- Equation 4.85 is just the signal model used to 
describe the received space-time coding MIMO SAR data in the elevation direction. 



4.5.3.3 Relations and Models for Azimuth Signal Processing 

In the azimuth direction, the space-time coding scheme results in equivalent phase 
centers when the radar moves from one position to another. This enables a coherent 
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Fig. 4.24 Extended DCPA 
scheme for the space-time 
coding MIMO SAR system 
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combination of the subsampled signals with the DPCA technique, which synthesizes 
multiple receive beams that are displaced in the along-track direction. It implies that 
we can broaden the azimuth beam from the diffraction-limited width, giving rise to 
an improved resolution without having to increase the system operating PRE 

Figure 4.24 shows the extended DCPA scheme for the space-time coding MIMO 
SAR system shown in Fig. 4.23, in which 6 transmit subarrays and 12 receive sub- 
arrays are employed. With the space-time coding scheme described previously, 
at azimuth time r/i the orthogonal signals X\(t) = [x a i(t), x a 2(t), x a ^(t)] T and 
Xb(0 = [xbi(t),Xb2(t),Xb3(t)] T are transmitted from the subarrays (1, 3, 5) 
and (2, 4, 6), respectively. Similarly, at azimuth time r\i the orthogonal signals 

X B« = [*w (0.*m (0,4, Ml 7 " and " X A<0 = ["**i(0, -** 2 (0, -*« 3 (0] r are 
transmitted from the subarrays (a, c, e) and (b, d,f), respectively. After beamforming 
processing in the elevation direction, it is possible to define the " two-way" phase 
centers as the mid-points between the equivalent transmit beam phase center and the 
equivalent receive beam phase center. All goes well as the radar samples are collected 
by the equivalent transmit and receive beams with phase centers collocated in the 
two-way phase centers. We aim to collect the radar returns at the time when their 
equivalent two-way phase centers occupy the at same spatial position, along with the 
platform trajectory. 

When the DCPA condition is matched, after beamforming processing in elevation 
we consider the geometry of equivalent beam phase centers, as shown in Fig. 4.25. 
The term a(r) is the equivalent beam instantaneous squint angle at azimuth time x, 
and R c is the closest range to a given point target when the MIMO SAR platform 
moves along its trajectory. For the signals transmitted at azimuth time x\ , the equiv- 
alent beam two-way phase centers between the equivalent transmit beam D and the 
equivalent receive beam A can be represented by 



Rda{t\) 



3Lg 

2 



-.2 



# c tan(a(Ti)) 



R* 



3L a 

2 



/? c tan(or(n)) 



(4.86) 
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Fig. 4.25 Geometry of 
equivalent beam phase 
centers 




for the signals transmitted in X , 
for the signals transmitted in X , 



where L a is the distance between subarrays in the azimuth direction. Similarly, for 
signals transmitted by the equivalent beam B at azimuth time r\ 2 and received by the 
equivalent beam C we have 



Rbc^i) 



La 



- tf e tan(o<(r 2 )) 



Rl 



La 



7? c tan(a(T2)) 



Kl 



(4.87) 
The corresponding equivalent phase difference between Rda( x i) an d ^Bc( r 2) is 



A*(Ti, T 2 ) = -^[RDA(Tl) - R BC (T2)] 

A 



(4.88) 



As the equivalent receive beam is assumed to be coincident in the far-field region 
and is of the same width as the equivalent transmit beam, for two successive pulse 
repeated intervals (PRI = 1/PRF). 



ot(Ti) «a(r 2 ) 



(4.89) 



As an example, assuming an X-band MIMO SAR system because this carrier 
frequency is popular in current spaceborne SAR systems with the following parame- 
ters: X — 0.03 m, R c — 700 km and L a — 4 m, Fig. 4.26 shows the corresponding 
equivalent phase difference A$(ti, r 2 ) as a function of the instantaneous squint 
angle a. We note that the phase difference is small and may be neglected during 
subsequent image formation processing. More importantly, in this case the channel 
responses H\ n and //2,„ at azimuth time ri and r 2 can be seen as being equivalent. 
Equation 4.82 can then be simplified into 



R\n' 
Rln 



(IS1I 2 
(IS1I 2 



|S 2 | 2 )ffi,„ 
|S 2 | 2 )#i,„ 



Nl,„ ~ N 2 ,„ 
Nl„ + N 2 ,n 



(4.90) 



Therefore, the signals received by the subantenna n can be recombined, so as to 
further improve its SNR performance. 
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Fig. 4.26 Equivalent phase 
difference AO(ti, T2) 
between two successive PRI 
transmission as a function of 
the instantaneous squint 
angle 
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The range from the equivalent MIMO SAR phase center to an arbitrary point 
target can then be represented by 



R{x) = JRl + v 2 (t - r ) 2 - 2v. 5 (r - r )fl cos(6>) 



(4.91) 



where Ro is the slant-range to broadside of antenna center, rjo is the time at target 
broadside, and 6 is the angle of the target off broadside. In SAR applications, cubic 
and higher order terms may be ignored in the azimuth phase history of targets. A 
general criterion for this condition to hold is that this term causes less than it/4 
excursion the phase over the aperture synthesis time. In this case, Eq. 4.91 can be 
simplified into 



R(r) = R - v, cos((9)(r - t ) + 



vjsin 2 (0)(T - r Q y 
2R 



(4.92) 



Therefore, for the case illustrated at the bottom of Fig. 4.24 the Doppler phase his- 
tories of the four received subarrays can be represented, respectively, by 



Vcm(t) 



;K™)+*(<-^)] 
^H*('-^) +s ('-^)] 



(4.93a) 
(4.93b) 

(4.93c) 

(4.93d) 
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(4.93e) 



— K R (-^) +s (-^)] 



(4.93f) 



(4.93g) 



(4.93h) 



where ir„ m {x) denotes the phase of the signal transmitted from the equivalent sub- 
array m and received by the equivalent subarray n. For example, Eq. 4.93c can be 
expanded into 



^Ac(r) = 



0.5L a 
2Rq — v s cos(#) I r tq 






2 ■ 2 

vt sin 



+ - 



(6) (e - 2=g» - ft,)' 



L/sin z (0) 
^1 



(4.94) 



In near-space vehicle-borne SAR applications, L a , 7?o and k are typically on the 
orders of 1 m, 20km and 10cm, respectively, hence the last term L a 2 sin 2 (6) / RqX is 
small and can be ignored. We then have 



2 / 0.5L a \ 
fAcir) = -R ( r - —^-\ = 

In a similar manner, we can get 



iAbs(t) 



itADir) = ifBc(r) = -R(t) 

A 



and 



<fcu(t) 



fBD(r) = -R 

A 



= 2 / 15M 

A- V v, / 



V04b(t) = iAba(t) =7^ 



(4.95) 

(4.96) 

(4.97a) 
(4.97b) 
(4.97c) 
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Fig. 4.28 MIMO SAR image formation processing via the rang-Doppler algorithm 



To construct one synthetic aperture, the system must operate with a PRF which 
leads, after combination of data streams, to a properly sampled synthetic aperture 
appropriate to the beamwidth of the system. Once the DPCA condition is met, 
the operating PRF of the system is always equal to the Nyquist sampling rate for 
the diffraction-limited beamwidth of the antennas. According to Eqs. 4.95-4.97, 
Fig. 4.27 illustrates the corresponding equivalent phase centers. In this way, for a 
fixed PRF value, the antenna length of the MIMO SAR can remain roughly con- 
stant with increasing N ea /2 (N ea is the number of equivalent beams in azimuth, 
and 2 results from the Alamouti coding scheme) the achievable azimuth resolution. 
Alternatively, for a fixed azimuth resolution, an increase in swath coverage can be 
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point target 




Fig. 4.29 Illustration of the range mitigation correction process 



obtained. If the DCPA condition is not matched, the gathered azimuth samples will 
be spaced nonuniformly. 

Figure 4.28 gives the range-Doppler (RD)-based image formation processing 
steps. As there are 6 orthogonal LFM signals transmitted by the transmit suban- 
tennas simultaneously, each receive subantenna will receive the total 6 signals. It is 
obvious that an echo with wide frequency width can be received by each subarray. 
As the transmitted signals are orthogonal, their matched filtering can be obtained 
separately with the corresponding reference functions. 

Thereafter, the combined data can be processed by a uniform image formation 
algorithm. First, range mitigation correction is employed, as shown in Fig. 4.29. 
Then, the focused image can be obtained through azimuth compression processing 
which is like normal monostatic SAR range-Doppler image formation algorithms. 



4.6 Conclusion 



Near-space vehicles provide a potential to high-resolution and wide-swath SAR 
imaging which is a contradiction for conventional SARs. This chapter describes 
the potential configurations including multichannel in azimuth, multiaperture in ele- 
vation, and space-time coding MIMO SAR. The corresponding system schemes, 
mathematical relations, and signal models are formed. It is anticipated that near- 
space vehicle-borne multiantenna SARs, especially MIMO SAR, enable SAR 
systems to be operated with high flexibility and reconfigurability, thus ensuring pre- 
viously unprecedented remote sensing performance. However, further research work 
is required. One future work is to develop high-precision image formation algorithms 
for MIMO SAR imaging. Another future work is synchronization processing. In this 
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chapter, we considered only monostatic MIMO SAR. Bi- and multi-static MIMO 
SARs are also feasible, but this configuration will result in technical challenges 
that lie in the synchronization between the transmitters and receivers, which include 
spatial synchronization, time synchronization, and phase synchronization. 



References 



1. Li, Z.F., Wang, H.Y., Su, T., Bao, Z.: Generation of wide-swath and high-resolution SAR 
images from multichannel small spaceborne SAR systems. IEEE. Geosci. Remote. Sens. Lett. 
2, 82-86 (2005) 

2. Currie, A., Brown, M.A.: Wide-swath SAR. IEE. Radar. Signal. Process. 139, 122-135 (1992) 

3. Curlander, J.C., McDonough, R.N.: Synthetic Aperture Radar: Systems and Signal Processing. 
John Wiley & Sons, Inc. (1991) 

4. The paradigm shift of effects-based space: near-space as a combat space effects enabler. 
http://www.airpower.au.af.mi. Accessed May 2010 

5. Suess M., Grafmuller B., Zahn R.: A novel high resolution, wide swath SAR system. In: 
Proceeingsof IEEE Int Geosci Remote Sens Symposium, Sydney, Australia 1013-1015 (2001) 

6. Suess M., Zubler M., Zahn R.: Performance investigation on the high resolution, wide swath 
SAR system. In: Proceedings of European Synthetic Aperture Radar Conference, Cologne, 
Germany 187-191 (2002) 

7. Heer C, Soualle E, Zahn R., Reber R.: Investigations on a new high resolution wide swath SAR 
concept. In: Proceedings of IEEE Int Geosci Remote Sens Symp, Toulouse, France 521-523 
(2003) 

8. Gebert N., Krieger G., Moreira A.: High resolution wide swath SAR imaging — system per- 
formance and influence of perturbations. In: Proceedings of Int Radar Symp, Berlin, Germany 
1-5 (2005) 

9. Gebert, N., Krieger, G.: Azimuth phase center adaptation on transmit for high-resolution wide- 
swath SAR imaging. IEEE. Geosci. Remote. Sens. Lett. 6, 782-786 (2009) 

10. Gebert N., Krieger G., Younis M., Bordoni E, Moreira A.: Ultra wide swath imaging with 
multi-channel ScanSAR. In: Proceedings of IEEE Int Geosci Remote Sens Symp, Boston, 
Massachusetts 21-24 (2008) 

11. Younis M., Bordoni E, Gebert N., Krieger G.: Smart multi-channel radar techniques for space- 
borne remote sensing. In: Proceedings of IEEE Int Geosci Remote Sens Symp, Boston, Massa- 
chusetts 278-281 (2008) 

12. Stiles J., Goodman N., SiChung L.: Performance and processing of SAR satellite clusters. In: 
Proceedings of IEEE Int Geosci Remote Sens Symp, Honolulu, Hawaii 883-885 (2000) 

13. Goodman, N., Lin, S., Rajakrishna, D., Stiles, J.: Processing of multiple-receiver spaceborne 
arrays for wide-area SAR. IEEE. Trans. Geosci. Remote. Sens. 40, 841-852 (2002) 

14. Aguttes J.P: The SAR train concept: An along track formation of SAR satellites for diluting 
the antenna area over N smaller satellites, while increasing performance by N. In: Proceedings 
of 55th Int Astronautical Congress, Vancouver, Canada 919-925 (2004) 

15. Aguttes J.P: The SAR train concept: required antenna area distributed over N smaller satel- 
lites, increase of performance by N. In: Proceedings of IEEE Int Geosci Remote Sens Symp, 
Toulouse, France 542-544 (2003) 

16. Griffiths H., Mancini P.: Ambiguity suppression in SARs using adaptive array techniques. In: 
Proceedings of IEEE Geoscience and Remote Sensing Symposium, Espoo, Finland 1015-1018 
(1991) 

17. Callaghan, G.D., Longstaff, I.D.: Wide-swath spaceborne SAR using a quad-element array. 
IEE. Radar. Sonar. Navig. 146, 159-165 (1999) 



96 4 Near-Space Vehicles in High-Resolution Wide-Swath Remote Sensing 

18. Fischer C, Heer C, Krieger G., Werninghaus R.: A high resolution wide swath SAR system. In: 
Proceedings of European Synthetic Aperture Rdar Conference, Dresden, Germany \-A (2006) 

19. Li, Z.F., Bao, Z., Wang, H., Liao, G.S.: Performance improvement for constellation SAR using 
signal processing techniques. IEEE. Trans. Aerosp. Electron. Syst. 42, 436-452 (2006) 

20. Li, Z.F., Bao, Z.: A novel approach for wide-swath and high-resolution SAR image generation 
from distributed small spaceborne SAR systems. Int. J. Remote. Sens. 27, 1015-1033 (2006) 

21. Jain, A.: Multibeam synthetic aperture radar for global oceanography. IEEE. Trans. Antenna. 
Propag. 27, 535-538 (1979) 

22. Jean, B.R., Rouse, J.W.: A multiple beam synthetic aperture radar design concept for geoscience 
applications. IEEE. Trans. Geosci. Remote. Sens. 21, 201-207 (1983) 

23. Goodman N., Rajakrishana D., Stiles J.: Wide swath, high resolution SAR using multiple 
receiver apertures. In: Proceedings of IEEE Int Geosci Remote Sens Symposium, Hamburg, 
Germany 1767-1769 (1999) 

24. Krieger, G., Moreira, A.: Spaceborne bi- and multistatic SAR: potential and challenges. IEE. 
Radar. Sonar. Navig. 153, 184-198 (2006) 

25. Wang WQ (2010) Bistatic Synthetic Aperture Radar Synchronization Processing. In: Kouemou 
G (ed) Radar Technology. In-Tech Press, India 

26. Wang, W.Q.: Multi-Antenna Synthetic Aperture Radar Imaging: Principles and Applications, 
In Chinese. National Defense Industry Press, Beijing (201 1) 

27. Krieger, G., Gebert, N., Moreira, A.: Multidimensional waveform encoding: a new digital 
beamforming technique for synthetic aperture radar remote sensing. IEEE. Trans. Geosci. 
Remote. Sens. 46, 31^16 (2008) 

28. Younis M., Venot Y., Wiesbeck J.: Digital beam forming on-receive-only for radar applications. 
In: Proceedings of German Radar Symposium, Bonn, Germany 213-217 (2002) 

29. Younis M.: Digital beam-forming for high resolution wide swath real and synthetic aperture 
rdar, Dissertation, Karlsruhe, Germany (2004) 

30. Krieger G., Fiedler H, Rodriguez-Cassola M., Hounam D., Moreira A.: Analysis of system 
concepts for bi-and multi-static SAR missions. In: Proceedings of IEEE Geosci Remote Sens 
Symposium, Toulouse, France 770-772 (2003) 

3 1 . Krieger, G., Gebert, N., Moreira, A. : Unambiguous SAR signal reconstruction from nonuniform 
displaced phase center sampling. IEEE. Geosci. Remote. Sens. Lett. 1, 260-264 (2004) 

32. Krieger G., Gebert N., Moreira A.: SAR signal reconstruction from non-uniform displaced 
phase center sampling. In: Proceedings of IEEE Int Geosci Remote Sens Symposium, Anchor- 
age, Alaska 1763-1766 (2004) 

33. Krieger G., Gebert N., Moreira A. (2004) Digital beamforming and non-uniform displaced 
phase centre sampling in bi-and multistatic SAR. In: Proc of European Synthetic Aperture 
Radar Conf, Ulm, Germany 563-566 

34. Gebert N., Krieger G., Moreira A.:SAR signal reconstruction from non-uniform displaced 
phase centre sampling in the presence of perturbations. In: Proceedings of IEEE Int Geosci 
Remote Sens Symposium, Seoul, Korea 1034-1037 (2005) 

35. Gebert N., Krieger G., Moreira A.: High resolution wide swath SAR imaging with digital 
beamforming-performance analysis, optimization and system design. In: Proceedings of Euro- 
pean Synthetic Aperture Radar Conference, Dresden, Germany 341-344 (2006) 

36. Gebert N., Krieger G., Moreira A.: Digital beamforming for HRWS-SAR imaging system 
design, performance and optimization strategies. In: Proceedings of IEEE Int Geosci Remote 
Sens Symposium, Denver, Colorado 1836-1839 (2006) 

37. Gebert N., Krieger G., Moreira A.: Multi-channel ScanSAR for high-resolution ultra- 
wide-swath imaging. In: Proceedings of European Synthetic Aperture Radar Conference, 
Friedrichshafen, Germany 79-82 (2008) 

38. Claassen J.P., Eckerman J.: A system concept for wide swath constant incident angle coverage. 
In: Proceedings of Synthetic Aperture Radar Technology Conference, Las Cruces, New Mexico 
41-59(1978) 



References 97 

39. Bellettini, A., Pinto, M.A.: Theoretical accuracy of synthetic aperture sonar micronavigation 
using a displaced phase-center antenna. IEEE. J. Oceanic. Enginneer. 27, 780-789 (2002) 

40. Lombardo, P., Colone, E, Pastina, D.: Monitoring and surveillance potentialities obtained by 
splitting the antenna of the COSMO-SkyMed SAR into multiple sub-apertures. IEE. Proc. 
Radar. Sonar. Navig. 153, 104-116 (2006) 

41. Wang, W.Q., Cai, J.Y., Peng, Q.C.: Conceptual design of near-space synthetic aperture radar 
for high-resolution and wide-swath imaging. Aerosp. Sci. Technol. 13, 340-347 (2009) 

42. Cumming, I.G., Wong, F.H.: Digital Processing of Synthetic Aperture Radar Data. Artech 
House, Boston (2005) 

43. Li, F.K., Johnson, W.T.K.: Ambiguities in spaceborne synthetic aperture radar systems. IEEE. 
Trans. Aerosp. Electron. Syst. 19, 389-397 (1983) 

44. Krieger G., Moreira A.: Potentials of digital beamforming in bi-and multistatic SAR. In: Pro- 
ceedings of IEEE Geosci Remote Sens Symposium, Toulouse, France 527-529 (2003) 

45. Younis, M., Fischer, C., Wiesbeck, W.: Digital beamforming in SAR systems. IEEE. Trans. 
Geosci. Remote. Sens. 41, 1735-1739 (2003) 

46. Gebert N., Krieger G., Moreira A.: Digital beamforming for HRWS-SAR imaging. In: Pro- 
ceedings of IEEE Geosci Remote Sens Symposium, Denver, USA 1836-1839 (2006) 

47. Gebert, N., Krieger, G., Moreira, A.: Digital beamforming on receive: techniques and optimiza- 
tion stragies for high-resolution wide-swath SAR imaging. IEEE. Trans. Aerosp. Electron. Syst. 
45, 564-592 (2009) 

48. Younis M., Patyuchenko A., Huber S., Bordoni E, Krieger G.: Performance comparison of 
reflector-and planar-antenna based digital beam-forming SAR. Int J Antenna Propag (2010) 
doi:10.1 155/2009/614931 

49. Huber S., Younis M., Patyuchenko A., Krieger G.: Digital beam forming techniques for space- 
borne reflector SAR systems. In: Proceedings of European Synthetic Aperture Radar Confer- 
ence, Aachen, Germany 962-965 (2010) 

50. Wang, W.Q., Peng, Q.C., Cai, J.Y: Waveform-diversity-based millimeter-wave UAV remote 
sensing. IEEE. Trans. Geosci. Remote. Sens. 47, 691-700 (2009) 

51. Wang, X.Q., Xiao, Q., Chen, Y.Q., Zhu, M.H.: The SNR study of the wide-swath SAR basing 
on elevation multi-receiver. J. Electron. Info. Technol. 29, 2101-2104 (2007) 

52. Krieger, G., Gebert, N., Moreira, A.: Multidimensional waveform encoding: a new digital 
beamforming technique for synthetic aperture radar remote sensing. IEEE. Trans. Geosci. 
Remote. Sens. 46, 31^15 (2008) 

53. Wang W.Q.: Applications of MIMO technique for aerospace remote sensing. In: Proceedings 
of IEEE Aerospace Conference, Big Sky, MT 1-10 (2007) 

54. Zhuge, X.D., Yarovoy, A.G.: A sparse aperture MIMO-SAR-based UWB imaging systems for 
concealed weapon detection. IEEE. Trans. Geosci. Remote. Sens. 49, 509-518 (2011) 

55. Cristallini, D., Pastina, D., Lombardo, P.: Exploiting MIMO SAR potentialities with efficient 
cross-track constellation configurations for improved range resolution. IEEE. Trans. Geosci. 
Remote. Sens. 49, 38-52 (2011) 

56. Wang W.Q.: Space-time coding MIMO-OFDM SAR for high-resolution imaging. IEEE. Trans. 
Geosci. Remote. Sens. (2011) doi:10.1109/TGRS.2011.2116030 

57. Picciolo M.S., Griesbach J.D., Gerlach K.: Adaptive LFM waveform diversity. In: Proceedings 
of IEEE Radar Conference, Rome, Italy, 1-6 (2008) 

58. Levanon, N., Mozeson, E.: Radar Signals. Wiley-IEEE Press, New York (2004) 

59. Kim J.H., Ossowska A., Wiesbeck W.: Investigation of MIMO SAR for interferometry. In: 
Proceedings of 4th European Radar Conference, Munich Germany 51-54 (2007) 



Chapter 5 

Near- Space Vehicles in Ground 

Moving Target Indication 



Abstract The last but not the least application of near-space vehicles in remote 
sensing is the Ground moving target indication (GMTI). In this chapter, we aim at 
MIMO SAR-based solution for GMTI including image formation processing. GMTI 
using MIMO SAR is interesting because the target indication performance can be 
improved by MIMO SARs' larger virtual aperture. The minimum detectable velocity 
of a target can be improved by both the larger virtual aperture size and the longer 
coherent processing interval used by MIMO SARs. 

Keywords Near-space • Ground moving target indication (GMTI) • MIMO SAR • 
Along-track interferometry (ATI) • Multi-antenna SAR • Fractional Fourier 
transform (FrFT) 



5.1 MIMO SAR with Multi- Antenna in Azimuth 

GMTI is of great important for surveillance and reconnaissance [1-6]. One repre- 
sentative GMTI technique is the along-track interferometry (ATI) SAR, which was 
initially proposed for estimating the radial velocity of ground moving targets, start- 
ing from interferometric phase measurements [7-9]. However, ATI SAR does not 
take into account the fact that the stationary clutter unavoidably corrupts the inter- 
ferometric phase of the target depending on its signal-to-clutter environment [10]. 
Consequently the imaged moving targets may be displaced in azimuth according to 
its radial velocity and superimposed upon clutter at a wrong location. Moreover, the 
estimated velocity may be ambiguous. To overcome this disadvantage, a method of 
moving target detection and location with three-frequency three-aperture ATI SAR 
was proposed in [1 1], but the parameter estimation and image formation algorithms 
are not discussed. Differently, here we discuss MIMO SAR-based GMTI approaches. 
The multiple antennas in MIMO SAR systems can be placed in elevation and/or 
azimuth, but we consider only multiple antennas in azimuth. Figure 5.1 shows the 
geometry of the MIMO SAR with multi-antenna in azimuth for GMTI applications. 
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Fig. 5.1 Geometry of 
MIMO SAR with 
multi-antenna in azimuth 




The SAR platform is moving along jt-axis at a constant speed of v s and at an altitude of 
h s . Suppose the transmit array has N, co-located antennas whereas the receive array 
has N r co-located antennas. The waveforms used in any two transmit antennas must 
be orthogonal, so that their echoes can be separated in subsequent signal processing. 
For each of the N t transmitted signals, there are N r coherent returns; hence, there 
will be Nt x N r different returns, each is weighed by different amplitude and phase. 
Each receive antenna signal is separately demodulated, digitized, and stored. Next, 
multiple pairs of equivalent transmitter-to-receiver phase centers can be formed by 
a jointly spatiotemporal processing. This provides a potential solution to GMTI 
processing. 

Also, we make the assumptions of far-field, flat earth, free-space, and single- 
polarization for our model. Suppose there are three along-track antennas. The range 
history from the central antenna to a given moving target located at {xq, yo) with 
a velocity of (v x , v y )(v x and v v denote the velocity elements in x-axis and y-axis, 
respectively) is represented by 



R 2 (x) = J[x + (y x - v,)t] 2 + (y + v v r) 2 + h 2 
[x Q + (v, - v.Jr] 2 + (vyt) 2 



RO + VvT + 



2*0 



where r is the azimuth slow time. The reference range Rq is 



R = Jy*+hl 



(5.1) 



(5.2) 



Similarly, the range histories from the target to the left antenna and right antenna 
are represented, respectively, by 



Rl(t) ^ RQ + Vyt + 



[xp + d a + (v s - v x )r] 2 + (v y T) 2 
2R 



(5.3) 
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„,. „ , [x -d a + (v s - V x )t] 2 + (Vyt) 2 

Ri(r) « Ro+Vyt + — (5.4) 

where d a is the separation distance between two neighboring antennas in azimuth. 
Without loss of generality, amplitude terms are ignored in the following discussions. 
Suppose s m (?) is the signal transmitted by the mm (m e (1,2, 3)) antenna with carrier 
frequency f c ,m, the echoes received by the nth (n e (1, 2, 3)) receive antenna can 
then be expressed as 



/ R m (r) + R n (r) \ 



\ C 

R m (r) + R n (r) 



x exp j2jtf ( 



co 



(5.5) 



where t and f c denote the range fast time and radar carrier frequency, respectively. As 
the transmitted waveforms are orthogonal, after separately range pulse compression, 
we have 

s n ,m(t, r) = exp (-;'27r/ Cjm T„ im ) • sinc[ir B r (t - r„, m )] (5.6) 

where x nm is the signal propagation time from the mth transmit antenna to the nth 
receive antenna 

R m (r) + R n (r) 

co 

The term s nm (t, r) denotes the matched filtering output at the nth receive antenna 
and matched to the mth transmit antenna waveform. This is just the signal model of 
the MIMO SAR with multi-antenna in azimuth, which means that there are a total 
of 3 x 3 matched filter outputs (or sub-images) resulting from a spatial convolution 
of the transmit and receive arrays. 

As the multiple antennas are displaced in azimuth direction, it is effective to define 
"two-way" phase center as the mid-point between the transmit and receive phase cen- 
ters. This provide a potential clutter suppression, like the DPCA technique. The clut- 
ter cancelation can be performed by subtracting the samples of the returns received 
by two-way phase centers located in the same spatial position, which are temporally 
displaced. The returns corresponding to stationary objects like the clutter coming 
from natural scene are cancelled, while the returns backscattered from moving tar- 
gets have a different phase center in the two acquisitions and remain uncanceled. 
Therefore, all static clutter are cancelled, leaving only moving targets. 



5.2 MIMO SAR-Based GMTI Processing 

Consider the three azimuth-antennas shown in Fig. 5.2. For the mth transmit antenna, 
the range pulse compression results at the three receive antennas are represented, 
respectively, by 
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Fig. 5.2 Scheme of MIMO 
SAR with three 
azimuth-antennas for GMTI 
applications 



Tx/Rx 



^ 




A 




si, m (t, r) = exp (-j2nf c , m T\ <m ) ■ sinc[jr B r (t - n, m )] (5.8a) 

sz,m(t, r) = exp (-j2nf c , m T2 <m ) ■ sinc[jrB r (f - r 2 , m )] (5.8b) 

*3,m0, r) = exp (-;'27r/ c , m T3, m ) ■ sinc[7rB r (f - t 3i ,„)]. (5.8c) 

To cancel stationary clutter, we perform DPCA processing 

SC m ,2l(t, t) =S2,m(t,t) -si, m (f,r+ At) ■ G (5.9) 

sc m ,3 2 (t, r) = S3, m (t, x - At) ■ G - S2, m (t, r) (5.10) 

where At = d a /(2v s ) is the relative azimuth time delay between two neighboring 
antennas, and Go is used to compensate the corresponding phase shift between two 
antennas 



Go = exp 



/_ 2^\ 
V J 4R Q k)' 



(5.11) 



Since there is 
A = sinc[7r5 r (f - n, m )] « sinc[jrfl r (f - T 2 , m )] » sinc|>rfl r (f - T 3 , m )]. (5.12) 
The Eqs. 5.9 and 5.10 can be derived into 



sc m> 2\(t, t) = A ■ exp (-;2jr/ c , m T 2 , m ) ■ exp 
1 — exp 



(-^' Ar ) 



'(-^ AT )] 



(5.13) 



•scm,32(* 



, t) = Ao-exp(-j27r/ c , m T 2 , m ) -exply— VyArj • 1 - exp f — /'— v y At j 

(5.14) 



5.2 MIMO SAR-Based GMTI Processing 
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where k m is the wavelength. It can be noticed that, if v y =0, there will bs\sc m> 2i(t, t) 
= 0| and\sc m T,2(t, t) = 0|; hence, the clutter is successfully canceled by this method. 
The next problem is to detect the moving targets. 

To estimate the moving target's radial velocity v y , we perform 



. . sc m ,32(?, r) 

sc m , 32-21 if, r) = = exp 

■sc m ,2iu, r) 



{ j4 i VyAX } 



We then have 



(An \ 

■^1,32-21 (?, t) = exp I j— v y Ar I 

(An \ 

sc 2 ,32-2iit,T) = exp I j — v v At I 

/ 4tt \ 

*c 3 ,32-2l(f>T) = expl ;'— v v Ar I. 



Next, we can get the following interferometry phases 

AOi = arg{ic 2 ,32-2i(f, t) ■ sc* 32 _ 2i (t,r + At)} 



exp[;4*v,Ar(l-l)] 

arg{5C 3 , 32-21 it, T) ■ SC2 i32 _21 < 

exp[^v y Ar(l-l)] 



A0 2 = arg{5c 3 , 32-21 (f, r) ■ sc^ 32 _2i (*, t + At)} 



A$ 3 =arg{jc 3i 32_2i(f, r) ■ sc] 1 32 _ 21 (/, r + 2Ar)} 
= exp[j4^Ar(i-i)]. 



(5.15) 

(5.16a) 
(5.16b) 
(5.16c) 



(5.17) 



(5.18) 



(5.19) 



The radial velocity v v can then be calculated from Eqs. 5.17 to 5.19. Taking 
Eq. 5.19 as an example, we can get 



V s A<$>3Xi\t, 



v. s A0 3 A.iA 3 



(5.20) 



AnvsAriki - X3) 2ndik\ — k$) 
The corresponding unambiguous speed estimation range v v „„, is determined by 

v s k-^k\ v s X^,X\ 



2d a ik\ - k 3 ) 



< v v < 



2d a (ki -k 3 )' 



(5.21) 



It can be noticed that the unambiguous speed estimation range has been signifi- 
cantly extended when compared to the conventional single-carrier detection method. 
This is particularly valuable for GMTI applications. However, it is necessary to be 
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noted that, to avoid possible blind speed problem in target detection, all the three 
interferometry phases should be used to calculate the radial speed. This means that we 
can get possible three speeds and the true speed is identified with jointly processing, 
so as to avoid loss detection of the moving targets. 

The double-interferometry method expressed in Eq. 5.20 has a good unambigu- 
ous speed estimate range which is an advantage for detecting high-speed targets, 
but its estimate precision may be unsatisfactory for some specific applications. To 
overcome this disadvantage, we firstly use the speed obtained from the previous 
double-interferometry method to de-ambiguity the targets. Next, after range mitiga- 
tion correction we can perform interferometry processing in range-Doppler domain 
for the same transmit antenna [11] 

A4>;„ = arg{sc m ,32(r, f d ) ■ < >21 (r, f d )} = 2 ™y< d ° , (5 . 2 2) 

We can then obtain the residual radial speed 

A m v,AcD' 



2jtd a 



(5.23) 



Therefore, both unambiguous speed estimate range and precision can be improved 
by this jointly processing algorithm. 

To estimate the moving targets' Doppler parameters, we transform the range 
equation expressed in Eq. 5.1 into 

_ , , _ „ *0 , >"0 V )' + X 0( V x - V,) l>y + (Vx ~ V s ) 2 2 

R 2 (r) « R + r-£- + - r + - — x\ (5.24) 

and substitute it into Eq. 5.13 which yields 

* C (T) = A\ expt-yjr^T 2 + 2f dc r) + </> ] (5.25) 



with 



Ai = A exp ( -j-^VyAtj ■ 1 - expl-j-^-VyAxj (5.26a) 

k d = -^—[(v s -v x ) 2 + v y ] (5.26b) 

2 
fdc = —-z-[yov x + x (v x - v s )l (5.26c) 

A2«0 

An / xi \ 

O = - T (* O+ ^L). (5.26d) 

Therefore, once the Doppler parameters k d and f dc are obtained with estimation 
algorithms (detailed in next section), the parameters of v x , yo , and Rq can be obtained 
by the Eqs. 5.2 and 5.26 because the h s in Eq. 5.2 is knowable from the inboard motion 
measurement or compensation sensors. 
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5.3 Simplified FrFT-Based Parameters Estimation 

To estimate the Doppler parameters k ( ] and fd c , we developed a simplified fractional 
Fourier transform (FrFT)-Based algorithm. 



5.3.1 Simplified FrFT Algorithm 

The conventional FrFT is a hybridized time-frequency transform algorithm. Its trans- 
form kernel is defined as [12] 



K a {t,IA) 



l = ^te e ;^_cote-;/«f M ca j a ^ kn 
S(f-n), a = 2kjt (5-27) 

8(t + fi), a = (2k+\)n 



where k denotes an integer, and a indicates the rotation angle in FrFT domain. 
This operation can be considered as a generalized form of Fourier transform that 
corresponds to a rotation over an arbitrary angle a — rir/2 with r e i)t. The forward 
and inverse FrFT of x(t) are defined, respectively, by 

00 

X«0*)= J x(f)K a (f,n)dt (5.28) 

— oo 

00 

x(t)= [ x a (iJ,)K- a (t,n)du. (5.29) 



The FrFT of a function x(t), with an angle a, can be computed in the following 
steps. 

Step 1. A product by a chirp: 

g(f)=.T(f)e-^' 2tan( 3 ) . (5.30) 

Step 2. A Fourier transform (with its argument scaled by csc(a)) or a convolution: 



g *(t) = g (f) © e -/'i' 2csc («) = J X (t)e- j ^- t)2csc(a) dt (5.31) 

-oo 

where © denotes a convolution operator. 
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Step 3. Another product by a chirp: 

f(t) = e - j ^ 2m( ^g*(t). (5.32) 

Step 4. A product by a complex amplitude factor: 



, /l — /cot(a)\ 
^[*(f)] = j( ^^ J/W- ( 533 ) 

As the Steps 3 and 4 are redundant for signal detection, we name the FrFT without 
the Steps 3 and 4 as the simplified FrFT (sFrFT) [13]. It is also a linear transform and 
continuous in the angle a. This paper just uses this sFrFT to estimate the Doppler 
parameters. 

5.3.2 Simplified FrFT-Based Estimation Algorithm 

To estimate the Doppler parameters, applying the sFrFT to the Eq. 5.25, we get 



X«(m) = ^oexpO'0i) 

Tp/2 

■ / exp 

-T p /2 



j2jif dc x + jirk d T- + j — cot(a) - jjix csc(a) + j'0i 



dr. 
(5.34) 



It arrives its maximum at [14] 



{a , Ao) = argmax|X„(/Li)| 2 (5.35) 



k d — -cot(ao) 

ao = ' Ar P--50) 

fdc = £ocsc(a ). 

This condition forms the basis of estimating the moving targets' parameters. In 
the sFrFT domain with a proper a, the spectra of any strong moving target will 
concentrate to a narrow impulse, and that of the clutter will be spread. If we can 
construct a narrow band-stop filter in the sFrFT domain whose center frequency 
around at the center of the narrowband spectrum of a strong moving target, then the 
signal component of this moving target can be extracted from the initial signal. With 
this method the strong moving targets can be extracted iteratively, thereafter the weak 
moving targets may be detectable. This method can be regarded as an extension of 
the CLEAN algorithm [15] to the sFrFT. 

Therefore, after canceling the stationary clutter, the identification of the moving 
targets can be implemented with sFrFT in the following steps: 
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K , fdc 



first interferometry processing 

V ■ ;C «,.32-21 



second interferometry processing 



K , fdc 



v . v . x, . y„ 



z 



moving targets image formation processing 



Fig. 5.3 The flow chart of the sFrFT and double-interferometry processing combined algorithm 



Step 1. Apply one sFrFT to the data in which the clutter has been canceled with 
different a, and from the maximal peak we can get the numerical estimation of 
{jX, a). 

Step 2. Apply F" to the same data, we have 



Xa(ti) = X&Qi)- 



(5.37) 



Step 3. After identifying the first moving target, we then construct a narrow band- 
stop filter M(p) to notch the narrow band-stop spectrum of this moving target 



X'Jfi) = X a (p)M(p). 



(5.38) 



Step 4. The filtered signals are then rotated back to time-domain by an inverse 
sFrFT. 

Step 5. Repeat the operations from Step 1 to Step 4 until all the desired moving 
targets are identified. 
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Table 5.1 Simulation 
parameters 



Parameters 



Values 



Carrier frequency 

Flying altitude 

Flying velocity 

Pulse duration 

Transmit signal bandwidth 

Antenna length of each aperture 

Position of the moving target 

Position of the stationary target 1 

Position of the stationary target 2 

Velocity of the target C 



1.25 GHz 

4000 m 

100 m/s 

5|ls 

100 MHz 

1 m 

(x = m, y = m) 

(x = 10m, y = 0m) 

(x =0mj = 10m) 

(y x = 10 m/s, v y = 5 m/s) 




-20 -10 10 20 
range direction [m] 



-10 10 20 30 
range direction [m] 



Fig. 5.4 Comparative processing results, a Before applying the algorithm, b After applying the 
algorithm 



Once the Doppler parameters of each target are obtained, we can get the v x . Now, 
the parameters (v* , v y ) and (x , y ) are all determined successfully. Next, the moving 
targets can then be focused with one uniform image formation algorithm, such as 
range-Doppler (RD) and chirp scaling (CS) algorithms. An example processing steps 
are given in Fig. 5.3. 



5.4 Simulation Results 



Unmanned vehicle-borne three-antenna stripmap MIMO SAR data from three point 
targets, one moving target and two stationary targets, are simulated using the parame- 
ters listed in Table 5. 1. Figure 5.4 shows the comparative processing results using the 
CS-based image formation processing algorithm. It can be noticed that the imaged 
moving target is overlapped with the two stationary targets. The moving target cannot 
be identified from this figure, because we cannot discern which is the moving one 
and which is the stationary one. Next, clutter cancellation is obtained by applying the 
proposed double-interferometry processing technique, leaving only moving targets 
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and a much simplified target detection problem. The sFrFT algorithm is then applied 
to estimate the Doppler parameters. Finally, the moving targets are focused by the 
image formation algorithm with the estimated Doppler parameters. 



5.5 Conclusion 

In this chapter, we proposed a scheme of MIMO SAR for ground moving targets 
detection and imaging. Although MIMO radar has received considerable attention 
in recent years, little work is related to MIMO SAR remote sensing. Our approach 
employs MIMO along-track antenna configuration and waveform diversity. After 
cancelling clutter and estimating the along-track velocity of the moving targets 
with the proposed double-interferometry processing algorithm and estimating the 
Doppler parameters of the moving targets with the presented sFrFT technique, GMTI 
is obtained by a jointly processing algorithm. Finally, the moving targets are focused 
with one uniform image formation algorithm. In this way, both target location and 
target velocity can be acquired, and high-resolution moving targets' SAR imagery 
can also be obtained. Simulation results show its validity. 

As the conventional DPCA SAR detection performance is noise limited and the 
conventional ATI SAR detection is clutter limited, the proposed approach combines 
their advantages and overcomes their disadvantages; hence, this approach is more 
effective and robust than the conventional DPCA SAR-based GMTI approaches. In 
particular, it is not dependent on a target's across-track velocity component or its 
Doppler shift, which is difficult to determine due to insufficient freedom degrees. 
Moreover, it also allows the estimation of target's true azimuth position directly 
from its measured position in the final SAR images. Therefore, the MIMO SAR and 
double-interferometry processing combined GMTI approach is elegant and effective 
in moving target identification. 
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Chapter 6 
Summary 



Abstract Near-space does indeed seem to be an efficient solution to future remote 
sensing applications. The use of cost-effective near-space platforms can lead to the 
solutions that previously thought to be out of reach for remote sensing customers. 
However, there are several technological challenges. In this chapter, we discuss real- 
istic near-space remote sensing issues and possible future work. 

Keywords Near-space • Near-space remote sensing • High-precision imaging • 
Waveform diversity • Three-dimensional imaging • Near-space vehicle 



6.1 Realistic Near-Space Remote Sensing Issues 

Weather could be a significant problem if near-space vehicles are not furnished 
with reliable sensors for on-site meteorological data with which vehicle controllers 
can predict turbulence, icing, and violent gusts. The experience with high-altitude 
tropospheric operations from around-the-world balloonist teams and weather teams 
could be collected and codified to aid computer predictions at higher altitudes. 

Another limitation is the difficulty in getting vehicles to altitudes higher than 
just the lower reaches of near-space. Increasing vehicle altitudes greatly increases 
vehicle balloon size, quickly making them too difficult to manage on the ground and 
during launch. A rule of thumb is that for every extra 30 km of altitude to reach, 
the vehicle needs to have twice the volume. Vehicles quickly become unmanageable 
for repetitive and hasty military use if we try to consistently reach high altitudes. 
This is perhaps why the Air Force Scientific Advisory Board recently considered 
near-space vehicles higher than 30 km "not viable" in the near future. However, 
this unfortunately leaves near-space vehicles in the altitude region that will probably 
make them more susceptible to enemy attack. 

Another possible problem may be the loss of tactical or operational surprise. 
First, low-tech adversaries, such as terrorists and third world countries, may not 
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have the capability to easily track or monitor the vehicles. Also, even radar-capable 
adversaries may have difficulty in tracking vehicles at longer distances due to their 
slow speed and small RCS. Secondly, near-space vehicles will be loiter for a long 
time, and possibly be launched with such repetition that their existence will not 
provide the enemy any valuable intelligence. 



6.2 Future Work 

There are key technological issues which are still under research and require further 
development. A brief overview of issues related with near-space remote sensing is 
given here. 



6.2.1 High-Precision Imaging Algorithm 

One future work is to develop high-precision image formation algorithms for 
near-space vehicle-borne azimuth-variant BiSAR and multi-antenna SAR systems. 
For near-space vehicle-borne azimuth-variant BiSAR it is a mandatory to know 
both the transmitter-to-receiver distance and target-to-receiver distance to properly 
focus its raw data, but they clearly depend on the target height. In this case, conven- 
tional imaging algorithms may be not suitable to accurately focus the collected data. 
Azimuth- variant BiSAR data can be focused in time-domain and frequency-domain 
or hybrid-domain. As high-computational cost is required for time-domain process- 
ing algorithms, the frequency-domain or hybrid-domain processing algorithms are 
of great interest. However, their processing efficiency depend on the spectrum model 
of the echo signal. One equivalent velocity-based two-dimensional spectrum model 
was proposed in Chap. 3, but only the response of a point target was simulated. 
We are aware that much further work is required to develop practical algorithms for 
distributed or multiple targets imaging. Also, near-space vehicle-borne multi-antenna 
SAR imaging algorithms requires much further work. Since multi-antenna SAR is 
in its infancy, there is a lack of efficient imaging algorithms for them. 



6.2.2 Waveform Diversity Design 

Continuous frequency coverage in the OFDM-LFM waveform is assumed in 
Chap. 4; however, if there are frequency overlaps or frequency separations in the 
waveform, unwanted side lobes will be generated in subsequent digital beamform- 
ing processing. Figure 6.1 shows the impacts of frequency overlap and separation 
on final matched filtering results. A frequency overlap of 0.5B and a frequency 
separation of 0.5B are assumed, respectively. It can be noticed that the frequency 
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Fig. 6.1 Impacts of frequency overlap and separation on final matched filtering results: a-c the 
case with frequency overlap, d-f the case with frequency separation 



overlap and separation interval between the sub-chirp signals play an important role 
on the sidelobe performance. When continuous frequency coverage is used, a good 
sidelobe performance can be obtained. Otherwise, the sidelobe performance will be 
degraded. To further suppress the sidelobes, some other efficient algorithms should 
be developed. 



6.2.3 Three-Dimensional Imaging 



Another future work is near-space vehicle-borne three-dimensional (3D) imaging. 
Conventional SAR obtains its two-dimensional (2D) high-resolution images along 
azimuth and elevation dimensions by projecting the 3D distributed targets onto 
the 2D plane. Consequently it usually suffers from geometric distortions, such 
as foreshortening and layover. The strong shadowing effects caused by building, 
hills, and valleys may result in the information loss of the explored area. In order 
to overcome these disadvantages, 3D SAR imaging has become an urgent need 
[1]. Current 3D SAR imaging techniques include interferometry SAR (InSAR), 
curved SAR (CSAR) [2-4], and linear array SAR (LASAR) [5]. InSAR can create 
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Fig. 6.2 Potential antenna configurations for near-space vehicle-borne SAR 3D imaging 



interferometric images to measure the scene height, but it has no ability to cre- 
ate height resolution cells. CSAR utilizes curve platform trajectory to achieve 3D 
imaging, but it is difficult to provide a precise curve trajectory. LASAR can pro- 
vide high-resolution 3D images by combining the real- and synthetic-aperture radar 
imaging technology. Two downward-looking LASAR systems, DRIVE [6-8] and 
ARTINO [9-11], are being developed at ONERA and FGAN-FHR, respectively. 
However, the antenna arrays used in these two systems are SIMO (single-input and 
multiple-output) configuration. It is necessary to extend it into MIMO configuration. 
Figure 6.2 shows two potential configurations for near-space vehicle-borne SAR 3D 
imaging. 
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